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I  BACKGROUND  AND  PERSPECTIVE 


The  Federal  Aviation  Administration's  Systems  Research  and  Develop¬ 
ment  Services  and  National  Aviation  Facilities  Experimental  Center  (NAFEC) 

i 

are  engaged  in  developing  and  applying  computer-aided  systematic  method¬ 
ologies  to  the  national  network  of  (1)  airways  and  (2)  ground  station 
VORs  and  VORTACs  used  to  create  airways  and  provide  other  services. ^ 
Stanford  Research  Institute  was  awarded  Contract  No.  DOT-FA71WA-2547  on 
19  February  1971, .to  support  the  FAA  in  this  effort.  The  work  under  the 
contract  was  to  be  done  in  close  cooperation  with  a  NAFEC  team. 

The  primary  objective  of  the  contracted  work  was  the  development  of 
the  following  two  analytical  tools: 

•  A  systematic  methodology  and  associated  computer  program(s)  for- 
the  design  of  a  national  network  of  airway  routes  for  the  optimal 
movement  of  high-altitude  (above  18,000  ft.)  area  navigation 
traffic  flow.  This  is  referred  to  as  Tasks  2  and  3. 

•  A  methodology  and  associated  computer  programs  for  the  synthesis 
of  a  VOR/VORTAC  network  to  provide  navigational  aids  for  the 
related  airway  routes.  This  is  referred  to  as  Task  4. 

Task  1  is  concerned  with  preparing  suitable  traffic  demand  models 
on  the  basis  of  present  and  forecast  traffic  data.  These  models  are  to 
be  used  'as  computer  input  data  for  development  and  testing  of  the  Hi-RNAV 
networks.  A  discussion  of  Task  1  has  been  included  as  Appendix  G  of  this 
part  of  the  report,  the  rest  of  which  is  concerned  with  SRI's  accomplish¬ 
ment  to  date  in  connection  with  Tasks  2  and  3. 

Various  models  and  techniques  to  design  a  national  VORTAC  grid . are 
described  in  Part  2  of  the  report.' 
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The  question  of  combining  the  route  network  design  and  the  VORTAC 
grid  design  is  discussed  in  Sections  IV  and  V,  and  Appendix  I  of  this 
volume,  A  similar  discussion  is  included  in  Part  2  of  the  report  for 
the  sake  of  making  both  parts  of  the  report  self-contained. 
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II  SUMMARY 


A .  General 

This  draft  report  describes  and  explains  the  research  done  by  SRI 
during  the  period  from  the  awarding  of  the  contract  until  the  end  of 
August  1972  to  develop  a  computer-aided  air  route  network  design  method¬ 
ology.  The  methodology  developed  so  far  does  not  yield  a  network  design 
directly.  Rather,  it  consists  mainly  of  tools--in  the  form  of  mathematical 
models  and  techniques--to  assess  various  network  attributes  and  to  im¬ 
plement  various  design  alternatives.  A  complete  methodology  giving  a 
direct  route  network  design  requires  acceptable,  realistic  network  evalu¬ 
ation  criteria  based  on  various  network  attributes  and  expressed  in  con¬ 
crete  analytical  terms.  It  is  hoped  that  the  results  obtained  so  far  and 
presented  in  this  report  will  provide  a  basis  for  discussions  and  con¬ 
sultations  with  various  experts  in  the  field  of  aviation,  after  which  it 
is  expected  that  a  complete  design  methodology  can  be  accomplished. 

B.  Air  Route  Network  Design  Tools  Developed  to  Date 

The  various  tools  and  techniques  that  have  been  developed  so  far  are 
briefly  described  below.  These  tools  and  techniques  can  eventually  be 
incorporated  in  a  complete  design  methodology. 

1.  Establishment  of  Route  Intersection  Points  and  Angles 

Mathematical  relationships  and  associated  computer  programs 
have  been. developed  by  which  the  intersection  points  of  various  routes 
or  route  sections  and  the  intersection  angles  can  be  established.  The 


latitude  and  longitude  of  the  start  and  end  point  of  each  route  or  route 
section  must  be  specified.  The  mathematical  relationships  are  based  on 
assuming  the  earth  to  be  a  sphere.  The  transformation  from  geocentric 
to  geodetic  coordinates  can  easily  be  accomplished  through  standard 
transformation  relationships. 

2.  Establishment  of  Various  Network  Attributes 

Mathematical  models  and  relationships  and  the  associated  com¬ 
puter  programs  have  been  developed  by  which  several  network  attributes 
can  be  established,  e.g., 

•  Expected  number  of  potential  conflicts  per  unit  time  at 
various  intersections  along  each  route. 

•  Expected  duration  of  the  potential  conflicts  at  various 
intersections. 

•  Capacity  at  various  intersections,  based  on  a  specified 
number  of  acceptable  potential  conflicts  and  on  conflict 
duration. 

•  Maximum  value  of  the  flow  product  at  crossing  routes,  with 
an  indication  of  the  altitude  level  at  which  this  value 
occurs. 

•  An  indication  if  the  existing  flow  product  at  any  inter¬ 
section  exceeds  a  prespecified  safety  margin. 

»  •  Flight  mileage. 

•  Route  mileage. 

•  Intersection  density  count. 

•  Existing  flows  along  various  routes  and  route  sections,  and 
an  indication  if  any  flows  exceed  a  prespecified  limit. 

The  flow  rate  of  aircraft  at  various  altitudes  between  various  communities 
must  be  specified.  Note  that  safety  and  delay  in  an  air  traffic  control 
(ATC)  system  are  implicitly  accounted  for  in  the  potential  conflict  count 
and  the  conflict  duration,  since  a  system  producing  fewer  potential 
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conflicts  or  shorter  potential  conflict  durations  is  expected  to  be  safer 
and  to  result  in  less  delay. 

3.  Merging  and  Bending  Routes 

Techniques  and  associated  computer  programs  have  been  developed 
to  merge  closely  located  routes  on  the  basis  of  ATC  convenience  and  prac¬ 
tical  considerations.  The  possibility  of  bending  a  route  at  suitable 
places  to  bypass  congested  or  restricted  areas  has  been  considered.  A 
convenient  systematic  technique  to  implement  this  possibility  has  not  yet 
been  developed. 

4.  Computer  and  Plot  Programs 

A  computer  program  has  been  developed  in  Fortran  IV  language  to 
implement  various  steps  of  the  design  methodology  and  is  described  in  de¬ 
tail  in  the  User's  Manual,  "Computer  Aided  Traf f ic/A irway/VOR(TAC )  Net¬ 
work  Program,"  by  John  J.  Penick  and  Kathryn  N.  Sapios  dated  October 
1972.  A  program  has  also  been  developed  to  plot  the  air  routes  based  on 
the  output  of  the  computer  program.  A  description  of  the  plot  program 
is  included  in  the  User’s  Manual.  The  program  is  compatible  with  the 
IBM  7090,  Sigma  5,  and  IBM  System  360. 


Ill  PROJECT  SCOPE,  GROUND  RULES,  AND  PROBLEM  FORMULATION 


A.  Project  Scope  and  Ground  Rules 


During  the  early  stages  of  the  project,  the  NAFEC  and  SRI  team  mem¬ 
bers  held  several  discussions  to  establish  the  project  scope  and  ground 
rules  so  as  to  limit  the  required  efforts  to  manageable  bounds.  The 
scope  and  the  various  ground  rules  relating  to  Tasks  2  and  3  that  were 
tentatively  agreed  upon  between  NAFEC  and  SRI  are  summarized  below.  For 
further  details,  refer  to  Appendix  A. 

•  Traffic  Data  Base 

The  IFR  peak  day  traffic  tapes  will  be  used  as  the  data  base  for 
traffic  demands. 

•  Interface  with  Terminal  Areas 

Departure  and  arrival  points  will  be  assumed  to  be  located  50  nmi 
from  the  center  of  the  community.  Climb  and  descent  will  be  re¬ 
garded  as  terminal  problems  and  will  not  be  handled  explicitly 
in  the  design  of  the  route  network. 

•  Restricted  Airspace 

Restricted  airspace  will  be  considered  in  the  design  of  the 
route  network. 

•  Route  Widths 

Route  widths  will  be  assumed  to  be  constant  for  the  design  of 
the  route  network,  at  least  during  the  first  phase  of  the  project. 

•  Altitude  Utilization  Rules 

The  hemispherical  rule  will  be  used  as  a  starting  point  from 
which  other  altitude  utilization  rules  will  be  investigated  for 
possible  inclusion  in  the  Task  3  methodology. 

•  Atmospheric  Conditions 

Atmospheric ■ conditions,  i.e. ,  weather,  wind,  and  also  jet  stream, 
will  not  be  considered  explicitly  in  the  initial  design  of  the 
route  network. 
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•  Form  of  RNAV 


All  aircraft  will  be  considered  to  be  equipped  with  the  p  -  9 
type  of  RNAV  equipment  with  slant-range  correction. 

•  VQRTAC  Coverage 

For  the  purpose  of  route  network  design,  it  will  be  assumed  that 
full  VORTAC  coverage  is  available  throughout . the  United  States. 

•  High-Altitude  Level 

Only  air  route  altitudes  at  or  above  18,000  ft  will  be  considered 

B .  Identification  of  Design  Variables 

A  route  is  defined  by  the  horizontal  projection  that  gives  the  geo¬ 
graphic  location  of  various  points  and  by  its  altitude.  Thus  the  design 
of  a  route  network  requires  the  design  of  the  horizontal  profile  and  the 
altitude  for  various  routes.  The  horizontal  profile  is  generated  by 
bending  or  merging  the  route  with  another  route  or  diverting  it  from 
some  route;  in  some  cases,  the  direct  great-circle  path  between  two  com¬ 
munities  may  prove  to  be  the  desired  horizontal  profile.  Furthermore, 
the  use  of  RNAV  equipment  makes  it  possible  to  utilize  parallel  routes. 
Thus  the  three  variables  to  be  considered  for  the  design  of  a  route 
network  are: 

•  Route  bending  (merging,  diverging)  in  a  horizontal  plane 

•  Use  of  parallel  routes 

•  Assignment  of  altitude. 

C.  An  Initial  Problem  Formulation 

In  the  past,  air  routes  have  been  designed  on  an  as-needed  basis. 

The  present  effort  may  be  the  first  to  treat  air  route  design  as  a  system 
design  problem.  The  formulation  of  the  problem  as  a  system  design  prob¬ 
lem  has  required  some  simplifying  assumptions,  since  otherwise  the  prob¬ 
lem  quickly  becomes  intractable. 


After  considering  several  factors  that  affect  route  network  design, 
the  NAFEC  and  SRI  team  members  decided  to  formulate  the  problem  tenta¬ 
tively  in  the  following  way.  Given  the  traffic  demand  between  various 
community  pairs  in  terms  of  expected  aircraft  flows  per  unit  time,  de¬ 
sign  a  route  network  (i.e.,  select  the  horizontal  profile,  altitude,  and 
number  of  parallel  routes)  between  respective  community  pairs  such  that: 

•  ATC  constraints  are  fulfilled. 

•  A  suitable  combination  of  the  following  factors  is  accomplished 

-  Expected  number  of  potential  conflicts 

-  Expected  conflict  duration 

-  Capacity  of  air  routes 

-  Capacity  of  intersections 
Additional  route  length 

-  Flight  mileage 

-  Number  and  density  of  route  intersection  points 

-  Assignment  of  an  altitude  not  necessarily  that  desired  by  the 
user. 

It  is  expected  that  factors  omitted  or  inadequately  treated  in  the 
initial  formulation  can  be  introduced  later,  when  a  concrete  methodology 
for  the  tentative  formulation  has  been  developed. 

The  above  problem  statement  is  essentially  qualitative  in  nature. 

It  was  necessary  to  develop  analytical  models  for  all  the  features  indi¬ 
cated  so  that  a  concrete  and  quantitative  basis  for  comparing  various 
route  networks  could  be  established.  These  analytical  models  are  dis¬ 
cussed  in  Section  IV.  A  preliminary  methodology  to  change  various  de¬ 
sign  variables  and  test  the  resulting  network  attributes  is  presented 


in  Section  V. 


IV  THE  FAMILY  OF  NETWORK  ATTRIBUTE  MODELS 

A.  ATC  Constraints 

From  the  point  of  view  of  users  flying  between  Communities  i  and  j, 
the  best  route  is  the  direct  great-circle  path  between  i  and  j,  since 
this  is  the  shortest  route.  However,  there  frequently  are  great-circle 
routes  from  Community  i  to  Communities  j,  k,  1,  m  ...  such  that  the 
angles  between  Routes  ij  and  ik,  or  ij  and  il,  and  so  forth,  are  very 
small.  For  example,  the  direct  great-circle  route  from  New  York  to  Denver 
has  an  initial  heading  of  278. 13°  and  the  route  from  New  York  to  San 
Francisco  has  an  initial  heading  of  281.61°.  The  difference  between 
these  two  headings  is  only  3.48°.  Under  such  circumstances,  ATC  con¬ 
trollers  cannot  treat  these  routes  as  two  distinct  airways  for  a  con¬ 
siderable  length,  since  the  lateral  separation  minimum  is  not  fulfilled. 
Such  routes  can  be  merged  for  a  suitable  distance  to  a  point  where  the 
turn-off  angle  to  the  destination  is  sufficiently  large  that  both  the 
pilot  and  the  ATC  controller  can  easily  accomplish  it  within  a  reasonably 
short  time.  An  example  of  the  merging  of  two  routes  is  shown  in  Figure  1, 
using  plane  trigonometric  relationships.  A  detailed  discussion  and  ap¬ 
plication  of  spherical  trigonometric  relationships  is  included  in 
Appendix  B.  When  the  value  of  9  (see  Figure  1)  is  specified,  e.g.,  10°  '  ' 

or  15°,  it  is  easy  to  calculate  the  length  of  the  merged  part,  OC  =  x, 
in  terms  of  the  original  length,  &  =  OB,  and  the  initial  angle,  a.  As  . 
shown  in  Appendix  B,  the  length,  x,  and  the  increment  in  route  length, 

A  =  (OC  +  CB  -  OB),  are  given  by 

x  =  X  (cos  <y  -  sin  <y  cot  9)  •  (1) 

i  11 
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FIGURE  1  MERGING  OF  A  ROUTE  ALONG  ANOTHER  ROUTE 
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Appendix  B  also  shows  that,  for  a  departing  angle  0  of  10°,  the  increase, 
in  the  length  of  the  route  being  merged  does  not  exceed  0. 4  percent  of 
the  direct  route  length  and  that,  for  a  departing  angle  of  15°,  the 
excess  length  is  less  than  0.9  percent. 

B.  Expected  Number  and  Duration  of  Potential  Conflicts 

One  of  the  most  significant  attributes  of  air  route  network  and 
ATC  system  is  safety.  Unfortunately  the  expression  safety  is  not  easy 
to  quantify.  However,  it  is  agreed  that  a  flight  is  likely  to  be  safer 
than  another  flight  if  it  encounters  fewer  potential  conflicts  or  if 
the  potential  conflicts  are  of  shorter  duration.  Extending  this  approach 
to>  the  entire  network,  it  seems  reasonable  to  state  that  a  network  which, 
for  a  given  traffic  flow,  produces  fewer  expected  potential  conflicts  or 
shorter  total  potential  conflict  duration  is  safer  than  a  network  pro¬ 
ducing  more — at  least,  in  an  average  sense.  Furthermore,  every  potential 
conflict  requires  control  action  from  ATC  controllers.  Thus  ATC  workload 
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is  significantly  dependent  on  the  number  and  type  of  the  potential  con¬ 
flicts  to  be  resolved.  Typically,  ATC  controller  action  consists  of  di¬ 
verting  one  of  two  potentially  conflicting  aircraft,  which  means  the 
introduction  of  delay.  It  therefore  appears  that  safety,  ATC  workload, 
and  delay  are  all  heavily  dependent  upon  potential  conflicts  and  their 
duration.  Consequently,  an  analytical  model  giving  the  expected  number 
of  potential  conflicts  and  their  duration  would  be  a  very  useful  tool  in 
evaluating  a  network. 

Mathematical  models  concerned  with  potential  conflicts  at  an  inter¬ 
section  as  a  function  of  intersection  angle,  traffic  flow^  separation 
minimum,  and  average  aircraft  speed  have  been  developed  and  are  discussed 
in  Appendix  C.  A  summary  of  these  models  is  given  below  for  convenience. 
All  the  models  refer  to  the  two  intersecting  routes,  AB  and  CD,  shown  in 
Figure  2. 


FIGURE  2  TWO  INTERSECTING  ROUTES  SHOWING  FLOWS  OF  AIRCRAFT  PER  HOUR 

‘V  V 

1.  Expected  Number  of  Potential  Conflicts  per  Hour 

Let  the  radar  separation  minimum  be  X  miles  and  the  average 
speed  of  aircraft  along  Routes  AB  and  CD  be  V'mph.  Then  the  expected 
number  of  potential  conflicts  per  hour,  denoted  by  E  ,  is  given  by 
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f  f 

1  2  o' 

E  =  2  — —  X  sec  —  per  hour  ,  (3) 

•  ,  C  V  2 

where  f  and  f  are  the  flow  rates  of  aircraft  per  hour  along  AB  and  CD. 

1  2 

I 

2.  Expected  Conflict  Duration  Time  per  Conflict 

The  expected  duration  of  a  potential  conflict,  t ,  at  an  inter¬ 
section  is  given  by 

n  X  o'  _  • 

T  =  —  —  cosec  —  hours  .  (4) 

4  V  2 

Note  that  T  does  not  depend  on  f^  or  f ^ .  It  indicates  the  expected 
duration  of  potential  conflict  at  the  intersection,  which  is  a  function 
of  only  a,  X,  and  V. 

3.  Conflict  Duration  Time  per  Hour  at  the  Intersection 

The  total  time  per  hour,  H,  during  which  various  aircraft  will 
remain  in  conflict  if  nothing  is  done  by  ATC,  is  given  by 

2 

X 

H  =  tt  —  f  f  cosec  <y  .  (5) 

V  1  2 

The  usefulness  of  this  number  lies  in  indicating  the  type  of  conflict 
and  the  relative  effort  required  by  the  controller  to  resolve  it. 


C.  Capacity  of  an  Air  Route 

1 .  Capacity  Based  on  Average  Speed  and  a  Safety  Factor 

Let  V  be  the  average  speed  of  an  aircraft  flying  along  a  route, 
and  let  y  be  the  longitudinal  separation  minimum.  Then  the  capacity  of 
the  air  route  can  be  defined  as 
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(1  +  k)y 


where  k  is positive  constant  greater  than  or  equal  to  zero.  This  con¬ 
stant  is  usually  referred  to  as  a  safety  factor.  Thus,  for  V  =  600  and 
y  =  20,  the  capacity  of  an  air  route  with  zero  safety  factor  would  be 
30  aircraft  per  hour.  The  capacity  with  a  safety  factor  k  of  5  would  be 
six  aircraft  per  hour,  and  so  forth. 


2 .  Capacity  Based  on  a  Certain  Allowable  Number  of  Overtakings 

Let  a  single  route  be  used  by  various  types  of  aircraft  having 

speeds  V  ,  V  ,  V  ...  mph,  V  being  greater  than  V  ‘  .  Let  f  be  the 
1*  2’  3  ’  i  i+1  i 

flow  per  hour  of  aircraft  having  speed  V  .  We  can  define  the  average 

i 

capacity  of  the  route  to  mean  the  sets  of  flows  f  ,  f  ,  f  . . . , ' f or 

1  2  3  ' 

which  the  number  of  overtakings  of  aircraft  per  unit  length  of  route  per 
unit  time  does  not  exceed  a  specified  value  G  . 

t  T 


As  shown  in  Appendix  D,  the  number  of  overtakings  G  is  given 
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T  V 


1 

It 
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f  (V  -  V  )  f  (V  -  V  )  f  (V  -  V  ) 

2  1  2.31  3.41  4 

+  -  +  -  +  .  . 

V  V  .  V 

2  3  4 


f  (V  -  V  )  f  (V  -  V  ) 

3  2  3  4  2  4 

_f_  -  +  ... 

V  V 


+  ■  ■  ■  > 


(7) 


where  JL  is  the  route  length. 

t 

For  any  specified  value  of  G^  and  &  (e.g.,  two  overtakings  per 
hour  per  100  nmi),  the  capacity  of  the  air  route  is  the  set  of  combina¬ 
tions  of  the  values  of  *f  ,  f^,  ...  that  satisfies  Eq.  (7): 
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3. 


A  Special  Case 


If  f  ,  f  ,  f  ....  are  known  as  fractions  of  the  total  flow  f , 
l'  2 9  3’  ’ 

then  Eq.  (7)  can  be  rearranged  in  a  more  convenient  form.  Let  - 

f  =  k  f  . 
i  i 


Then, 


G  =  if 
T 


k 

1 

r k  (v  -  v  ) 
2  1  2 

k  (V  -  V  )  1 

3  13 

i  i 

k 

2 

rk  (V  -  V  )  1 

3  2  3 
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V 
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V 

L  2 

f  “T 

V 

3 

+  — 

V 

2 

-r 

V 

L  3  -1 

T 

,  (8) 


G  /i 

T 

f  =  .  /  — ■ - — ‘  aircraft  per  hour  , 

0(k.,V.) 

l  l 


(9) 


where  0(k  ,  V  )  is  equal  to  the  expression  shown  in  "braces"  in  Eq.  (8). 
i  i 

Equation  (9)  gives  the  capacity  of  a  single  air  route  in  terms  of  (G^/i), 

the  allowable  average  number  of  overtakings  per  hour  per  unit  length  of 

air  route.  The  speeds  of  various  aircraft  and  the  fraction  of  the  total 

flights  for  each  speed  are  assumed  to  be  known.  (Note:  It  is  assumed 

that  V  >  V  without  any  loss  of  generality.) 
i  i+1 

4.  Establishment  of  Parallel  Routes 

Parallel  routes  will  be  desirable  whenever  the  flow  along  a 

route  or  a  route  section  exceeds  a  prespecified,  limit  based  on  route 

capacity  considerations.  I f  we  let  R  indicate  the  allowable  maximum 

c 

flow  of  aircraft  per  hour  along  a  route  or  a  route  section,  then  the 

need  for  and  the  number  of  parallel  routes  can  be  established  by  dividing 

the  actual  flow  along  a  route  by  R  .  For  example,  let  the  flow  along  a 

c 

route  or  route  section  i  be  a  at  altitude  L.  Then 

•  ,  '  iL 

•  If  a  /R  is  less  than  1,  a  single  route  is  satisfactory 

'  iL  c 
at  level  L. 

•  If  a  /R  is  greater  than  1  but  less  than  2,-  two  parallel 

iL  c  ' 

routes  are  desirable  at  level  L. 
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Similarly  if  a  /R  is  greater  than  2  but  less  than  3,  three  parallel 

-  ..  :  *  iL  c 

routes  will  be  needed,  and  so  on. 

The  spacing  and  the  width  of  the  routes  will  depend  upon  the 
quality  of  the  available  VORTAC  coverage  and  accuracy  of  the  navigation 
equipment.  Suppose  an  air  route  width  W  nmi  is  established  as  a  nominal 
safe  route  width  in  the  region  where  the  need  for  parallel  route  has  been 
indicated.  The  center  lines  of  the  parallel  routes  can  then  be  laid 
(W+q)  nmi  apart  with  q  nmi  as  a  safety  zone.  The  route  width  is  dis¬ 
cussed  further  later  in  this  section. 

D.  Capacity  of  an  Intersection 

The  expected  capacity  at  the  intersection  of  two  routes  could  be 

defined  as  the  maximum  allowable  value  of  the  product  f  f  (where  f  and 

12  1 

f^  are  the  flows  of  aircraft  per  hour  on  the  two  routes)  under  one  of  the 
following  two  restrictions: 

(1)  The  expected  number  of  potential  conflicts  per  hour  does  not 
exceed  a  certain  specified  number. 

(2)  The  expected  time  that  the  aircraft  along  the  two  routes  vio¬ 
late  the  separation  minimum  does  not  exceed  a  specified 
fraction  of  an  hour. 


17 


E  V 
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(11) 


2  X  sec  — 

Equation  (11)  gives  the  capacity  at  an  intersection  as  a  func¬ 
tion  of  V,  X,  O' ,  and  some  specified  number  of  allowable  potential  con¬ 
flicts  E  per  unit  time, 
c 

2 .  Intersection  Capacity  Model  under  Restriction  2 
From  Appendix  C,  we  know  that 
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Equation  (13)  gives  the  capacity  in  terms  of  V,  X,  0/ ,  and  some 

specified  allowable  conflict  duration  time  H  per  hour.  For  further  de- 

i 

tails  of  these  mathematical  models,  the  reader  is  referred  to  Appendix  D. 

E .  Additional  Route  Length  per  Route  and  Route  Width 

In  the  present  air  route  system,  some  routes  between  community  pairs 
are  6  to  8  percent  longer  than  the  direct  great-circle  route  lengths. 

One  of  the  reasons  for  such  additional  lengths  is  that  the  present  routes 
are  based  on  travel  from  or  toward  a  VOR  or  VORTAC  facility,  and  the  lo¬ 
cation  of  VORTACs  is  frequently  not  in  line  with  the  direct  route.  In 
the^case  of  RNAV  routes,  this  restriction  is  not  present;  however,  even 
in  this  case,  some  routes  are  merged  or  bent  for  ATC  convenience  and  net¬ 
work  simplicity.  A  useful  attribute  for  route  network  design  would  be 
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the  length  added  to  each  route  because  of  merging  or  bending,  obviously 
a  simple  calculation.  The  extra  length  should  be  kept  to  a  minimum, 
since  users  would  prefer  as  direct  a  route  as  possible. 

Another  useful  attribute  to  compare  various  designs  of  a  route  net¬ 
work  (and  a  VORTAC  grid)  would  be  the  quality  of  coverage  provided  by 
the  VORTACs  to  various  routes  of  the  network.  A  convenient  method  of 
estimating  the  quality  of  coverage  is  to  compute  the  route  widths  along 
various  routes  using  the  route  width  diagram  in  FAA  Handbook  7110.18, 

Air  Traffic  Control  Service  for  Area  Navigation  Equipped  Aircraft,  Ap¬ 
pendix  D,  page  13  (27  February  1970).  Knowing  the  latitude  and  longitude 
of  the  VORTACs  supporting  various  routes  or  route  sections,  an  estimate 
can  be  made  of  the  route  widths  at  various  portions  of  the  route.  The 
percentage  of  total  route  length  supportable  with  a  prespecified  route 
width,  e.g.,  8  nmi,  could  then  be  used  as  a  measure  of  goodness  of  the 
design  of  the  route  network  and  the  VORTAC  grid. 

In  Appendix  I,  some  simplified  relationships  have  been  developed  to 
show  the  percentage  of  a  route  having  a  route  width  W  as  a  function  of 
the  orientation  of  the  route  with  respect  to  a  triangular  VORTAC  grid  of 
a  given  side  length  S.  The  percentages  are  calculated  for  each  route 
through  a  simple  subprogram. 

F.  Total  Flight  Mileage 

Flight  mileage  can  be  defined  as  the  product  of  the  length  of  a 
route  in  nmi  and  the  number  of  flights  per  unit  time  initiated  on  this 
route.  For  network  design,  it  is  not  sufficient  to  consider  just  route 
length;  it  is  also  necessary  to  take  into  account  how  frequently  each 
route  is  used.  An  infrequently  used  route  might  be  designed  with  a 
greater  percentage  of  extra  length  than  a  frequently  used  one.  This  as¬ 
pect  of  a  route  network  can  be  expressed  in  terms  of  the  total  flight 
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mileage  associated  with  the  route  network  for  some  given  traffic  demand, 
and  would  be  a  significant  network  attribute. 

The  calculation  of  total  network  flight  mileage  is  again  a  simple 
operation;  it  could  easily  be  implemented  through  a  computer  program. 

The  percentages  of  various  routes  having  a  route  width  W,  as  discussed 
above,  could  be  multiplied  by  the  aircraft  flows  along  the  respective 
routes;  the  percentage  of  flight  mileage  flyable  with  route  width  W 
could  also  be  computed  easily  thereby  providing  a  further  measure  of 
goodness  of  the  design. 

G.  Congestion  of  Route  Intersection  Points 

Congestion  of  several  route  intersection  points  in  a  small  area  is 
an  undesirable  feature,  from  the  points  of  view  of  both  the  ATC  control¬ 
lers  and  the  users,  e.g.,  pilots.  For  the  ATC  controllers,  it  means  in¬ 
creased  surveillance  and  control  efforts,  since  they  have. to  observe 
several  possible  locations  for  potential  conflicts.  For  the  users,  it 
means  frequent  maneuverings  to  avoid  conflicts  and  congestion.  A  measure 
that  indicates  the  location  per  unit  area  of  route  intersection  points 
would  be  a  useful  network  attribute.  One  method  of  evaluating  this  at¬ 
tribute  would  be  to  establish  the  maximum  number  of  route  intersection 
points  in  a  unit  area,  e.g.,  a  square  area  of  100-nmi  side.  This  could 
be  done  either  by  observation,  geometrically,  or  by  searching  the  network 
in  small  steps  in  a  systematic  manner,  using  a  computer  program.  A  sub¬ 
program  to  test  the  total  number  of  intersections  as  well  as  their  density 
in  a  square  region  of  specified  side  length  is  included  in  the  computer  . 
program. 
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V  THE  DESIGN  METHODOLOGY 

A .  The  Concepts  of  Threads,  Rope  Sections,  and  a 

Primary  Spider  Web  Network 

In  this  section,  we  present  some  concepts  agreed  upon  between  the 
NAFEC  and  SRI  team  members  during  various  working  sessions  in  the  early 
stages  of  the  project.  It  was  decided  that  during  the  experimental  de¬ 
velopment  of  a  route  network  design  methodology,  the  terms  "route"  or 
"airways"  should  be  avoided,  except  when  these  terms  are  used  as  defined 
in  current  FAR.  The  following  concepts  and  definitions  pertaining  to 
routes  and  airways  were  introduced  for  clarity  and  convenience. 

1.  Thread 

A  thread  is  the  shortest  distance  between  two  points.  These 
points  can  be  either  the  community  pairs  themselves  or  some  preselected 
points  in  the  vicinity  of  the  corresponding  community  centers.  A  de¬ 
scription  of  how  the  start  and  end  points  of  a  thread  could  be  established 

* 

is  presented  in  NAFEC Ts  Technical  Note  No.  8  (May  1971). 

2.  Primary  Thread  Direction 

A  primary  thread  direction  is  either  the  direction  of  one  of 
the  threads  connecting  two  community  centers,  or  a  fictitious  direction 
from  a  community  center  (or  an  associated  preselected  point)  along  which 


Bibliographies,  including  NAFEC  Technical  Notes,  are  given  at  the  end 
of  each  part  of  the  report. 
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one  or  more  adjacent  threads  will  be  made  to  lie  for  some  distance  for 
any  of  the  following  purposes: 

•  Simplifying  the  network 

•  Reducing  the  number  of  intersecting  routes 

•  Improving  safety  (e.g.,  reducing  the  number  of  potential 
conflicts) . 

3.  Rope  Section 

Any  section  of  a  route  that  consists  of  the  combination  of  two 
or  more  threads  is  termed  a  rope  section. 

4.  Primary  Spider  Web  Network 

The  network  resulting  from  connecting  by  threads  all  the  com¬ 
munity  pairs  exchanging  traffic  is  referred  to  as  the  primary  spider  web 
network.  Any  route  merging  or  bending  will  be  applied  to  this  primary 
spider  web  network.  Note  that  in  a  spider  web  network,  each  thread  be¬ 
tween  two  communities  is  a  great-circle  path. 


B .  Techniques  for  Merging  Threads  to  Form  Rope  Sections 

The  following  ground  rules  and  criteria  were  tentatively  agreed  upon 
between  the  NAFEC  and  SRI  team  members  for  the  establishment  of  route 
networks,  simplified  on  the  basis  of  some  physical  and  engineering 
considerations. 

Rule  1:  Where  applicable,  establish  a  few  primary  thread  directions 
from  a  reference  community  center,  such  that  the  angle  be¬ 
tween  two  adjacent  primary  thread  directions  is  at  least 
18°.  The  motivation  for  choosing  18°  is  that,  typically, 

50  nmi  of  travel  from  the  reference  community  is  needed 
for  jet  aircraft  to  obtain  an  altitude  above  18,000  ft. 

At  50  nmi  from  the  community  center,  the  minimum  angle 
between  two  threads  emanating  from  the  same  community  such 
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that  the  typical  airways  of  8  nmi  width  associated  with 
the  threads  will  cease  to  overlap  is  about  9°.  Thus 
threads  that  are  inclined  at  ±9°  or  less  with  reference  to 
a  primary  thread  can  be  treated  as  a  part  of  the  primary 
thread  for  the  initial  50  nmi.  Two  primary  threads  that 
are  inclined  at  least  18°  or  more  from  each  other  will  be 
mutually  exclusive  and  will  be  independent  of  each  other. 

Rule  2:  After  establishing  a  primary  thread  (or  its  direction) , 
merge  the  threads  in  its  vicinity  by  using  one  of  the 
following  approaches: 

•  Merge  to  meet  ATC  constraints  (i.e.,  minimum  departing 
angle  rule) . 

•  Merge  by  allowing  a  certain  percentage  mileage  increase 
for  the  threads  to  be  merged. 

Rule  3:  If  both  ends  of  a  thread  are  mergeable  along  different 
primary  threads,  use  one  of  the  following  approaches: 

•  Use  the  full  length  of  the  thread  independently  at 
each  end,  giving  two  independent  one-way  routes. 

•  Use  half  the  length  of  the  thread  for  one  end  and  the 
other  half  for  the  other  end. 


C .  Route  Merging  Based  on  Minimum  Departing  Angle 

Referring  to  Rule  1  above,  it  was  decided  that,  for  the  initial 
conception,  threads  will  be  merged  on  the  basis  of  the  minimum  departing 
angle  rule,  since  this  rule  is  relevant  to  ATC  operations  and  the  re¬ 
sulting  additional  length  is  not  excessive.  Consider  now  a  Community  i 
from  which  Threads  ij,  ik,  and  il  emanate.  Any  one  of  these  could  be 
selected  as  one  of  the  primary  thread  directions.  Suppose  Thread  ij  has 
a  heading  of  Y  .  and  that  it  is  selected  as  a  primary  direction;  then,  by 
Rule  1,  the  other  primary  directions  should  have  headings  of  y  .  +  A, 

7^  +  2A,  and  Yj  “  A,  y.  “  2A,  .  ..,  and  so  forth.  A  is  the  heading 

increment,  typically  15-20  degrees.  Each  of  the  threads  emanating  from 
Community  i  can  be  merged  partially  along  one  of  the  primary  headings 
Y  Y  .  +  A,  . , . ,  which  are  closest  to  the  thread.  The  same  technique 
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can  be  applied  to  all  other  communities.  By  selecting  various  combina¬ 
tions  of  reference  headings,  several  alternative  configurations  can  be 
studied  and  the  ones  that  appear  to  have  the  most  desirable  combination 
of  network  attributes  can  be  retained  for  detailed  study. 

In  this  connection,  it  should  be  noted  that  the  heading  increment 
A  and  the  departing  angle  9  should  be  so  selected  that  0  >  A/2,  other¬ 
wise  the  routes  making  an  angle  greater  than  0  but  less  than  A/2  with 
a  standard  heading  will  remain  unmerged.  For  example,  if  A  is  chosen 
as  20°  and  9  as  8  degrees,  a  route  making  an  angle  of  10°  with  a  standard 
heading  will  remain  unmerged. 

D.  Assigning  Altitudes 

It  was  established  during  the  early  stages  of  the  project  that  as¬ 
signment  of  altitude  is  not  a  practicable  design  variable.  As  shown  in 
•Appendix  E,  the  extra  cost  of  fuel  is  excessive  when  a  user  is  assigned 
an  altitude  other  than  the  one  he  has  chosen  on  the  basis  of  optimizing 
fuel  consumption  and  aircraft  performance.  Therefore,  altitude  assign¬ 
ment  as  a  design  variable  was  discarded.  It  is  assumed  that  the  user 
will  be  assigned  whatever  altitude  he  desires  within  the  framework  of 
the  hemispherical  rule  unless,  of  course,  other  considerations  apply, 
e.g.,  severe  weather  conditions.  If  the  hemispherical  rule  is  disre¬ 
garded,  possibilities  of  two-way  air  routes  having  the  same  altitude 
could  be  explored.  It  is  proposed  to  study  this  aspect  in  a  later  phase 
of  the  project.  The  only  practicable  design  variables  for  the  present 
preliminary  approach  are  thus  route  merging  and  bending,  and  route 
paral lei ing. 
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E .  Summary  of  the  Design  Methodology 


Below ,  we  summarize  the  design  methodology  accomplished  so  far,  for 
which  a  computer  program  has  been  developed  and  is  ready  in  an  operating 
form.  The  methodology  is  based  on  the  following  two  policy  decisions: 

1.  Routes  will  be  merged  using  the  minimum  departing  angle  rule. 

2.  Eastward  and  westward  routes  will  be  considered  separately. 

The  first  policy  decision  has  already  been  discussed.  The  second  policy 
decision  is  based  on  the  following  considerations: 

•  Since  the  Area  Navigation  techniques  offer  the  possibility  of 
establishing  independent  east  and  west  routes  (it  not  being 
necessary  to  fly  toward  or.  away  from  specific  VORTACS)  and  since 
the  altitude  bands  of  east  and  west  routes  are  mutually  exclusive, 
there  seems  to  be  no  reason  at  this  stage  to  mix  the  eastbound 
and  westbound  traffic  and  routes.  It  seems  proper  to  optimize 
the  east  and  west  routes  independently  as  a  first. step  and  then 
study  other  possibilities,  i.e.,  using  the  same  route  section  for 
both  eastbound  and  westbound  flow, 

•  Mixing  of  east  and  west  routes  would  require  use  of  more  than  one 
bend  in  the  route  when  it  is  merged  at  both  ends  along  different 
headings,  and  therefore  would  require  a  more  complex  programming 
effort.,  More  than  one  bend  in  a  route  may  not  be  desirable  from 
other  considerations  also. 

In  view  of  the  above,  it  was  decided  to  consider  the  eastward  and  west¬ 
ward  route  networks  separately  in  the  preliminary  development  of  the  pro¬ 
gram.  It  should  be  noted  that  with  suitable  modifications,  the  program 
can  be  used  to  consider  the  merging  of  half  route  lengths  along  respec¬ 
tive  headings  at  each  end.  We  now  describe  the  various  steps  of  the 
methodology: 

1.  Start  with  direct  great-circle  .threads  between  each  community 
pair: 

-  Each  thread  will  start  at  a  point  D-nmi  (the  assumed  terminal 
area  radius)  away  from  one  community  and  terminate  at  a  point 
D-nmi  away  from  another. 
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-  Given  the  latitude  and  longitude  of  two  points,  the  equation 
of  a  great-circle  path  is  a  standard  algebraic  expression: 

f(0>'M1>’l'1>92>+2>  =  0 

Establish  the  intersection  points  between  various  routes: 

-  A  test  whether  two  routes  intersect  is  available. 

The  latitude  and  longitude  of  each  intersection  point  can 
be  calculated,  as  well  as  the  angle  of  intersection. 

Calculate  the  network  attributes  of  interest: 

The  expected  number  of  potential  conflicts  and  their  dura¬ 
tion  along  each  route  per  unit  time. 

-  The  total  number  of  expected  potential  conflicts  and  their 
duration. 

The  number  of  route  intersection  points  in  the  horizontal 
plane  and  their  density  per  unit  area. 

-  The  total  flight  mileage. 

-  The  total  extra  flight  mileage  compared  with  the  mileage 
for  direct  routes. 

The  maximum  deviation  in  route  length  compared  with  the 
direct  length. 

Flow  along  various  routes  and  route  sections. 

Flow  product  at  various  intersections  and  intersection 
capacities. 

Select  a  reference  route  (or  a  reference  heading)  for  each  com¬ 
munity  and  establish  other  standard  headings  from  each  community, 
using  Rule  1  of  Section  V-B  and  elaborated  in  Section  V-C. 

Merge  the  threads  emanating  from  various  communities,  using 
the  minimum  departing  angle  rule,  and  compute  various  network 
attributes. 

Repeat  Step  4  by  selecting  other  suitable  reference  headings 
for  various  communities. 

Select  those  configurations  that  appear  to  have  the  most  de¬ 
sirable  combination  of  network  attributes  for  detailed 
evaluation. 
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F.  '  Applying  the  Methodology  to  Design  the  National  Network 

The  methodology  described  above  can  theoretically  be  applied  simul¬ 
taneously  to  any  number  of  communities  and  community  pairs.  However, 
the  computer  storage  requirements  and  computation  times  would  become 
prohibitive  if  excessively  large  numbers  of  community  pairs  are  con¬ 
sidered  at  the  same  time.  Also,  the  human  designer  would  find  it  ex¬ 
tremely  difficult  to  perceive  and  analyze  the  network  attributes  for  a 
very  large  number  of  routes,  intersections,  conflicts.,  etc.  Given  both 
these  human  and  computer  limitations,  some  form  of  decomposition  of  the 
design  process  is  essential. 

Fortunately,  the  nature  of  traffic  and  community  pair  combinations 
in  the  United  States  appears  to  have  some  inherent  features  that  can  be 
exploited  to  decompose  the  design  process  into  manageable  parts,  which 
can  be  completed  separately  and  then  combined  later  to  produce  a  national 
network.  We  start  with  the  following  observations  based  on  the  study  of 
the  peak  IFR  data  supplied  to  us  by  NAFEC. 

{ 

•  Total  number  of  domestic  flights  per  day  above  18,000  ft  is  on 
the  order  of  11,000. 

•  The  approximate  number  of  community  pairs  (threads)  exchanging 
various  levels  of  traffic,  and  the  total  number  of  flights  bn 
these  threads  are  shown  in  Table  1.  The  flights  are  also  ex-, 
pressed  as  a  percentage  of  total  flights. 

Table  1  indicates  that  a  relatively  small  number  of  community  pairs 
accounts  for  a  major  portion  of  the  traffic.  For  example,  if  only  the 
pairs  exchanging  10  or  more  flights  are  considered,  there  will  be  250 
pairs  (threads),  i.e.  about  11  percent  of  the  total  number  of  threads, 
and  these  serve  more  than  50  percent  of  the  total  traffic.  Similarly, 

500  community  pairs,  i.e,,  pairs  exchanging  5  or  more  flights  (23  percent 
of  the  threads),  account  for  over  70  percent  of  the  traffic.  Conversely, 
with  an  increase  o.f  only  17  percent  in  traf f ic--f rom  72  to  89  percent-- 
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Table  1 


/ 

NUMBER  OF  COMMUNITY  PAIRS,  TRAFFIC  LEVELS, 
AND  TOTAL  FLIGHTS  ABOVE  18,000  FEET 


Traffic  Level 

(exchange  of  flights 
per  day) 

Cumulative 

Number  of 

Community  Pairs 
(approximate) 

Sum  of 

Total  Flights 
( approximate) 

Cumulative 

Percentage  of 
Total  Flights 

20  or  more 

100 

4,110 

37 

15  or  more 

150 

5,010 

46 

10  or  more 

250 

6,210 

57 

5  or  more 

500 

7,960 

72 

2  or  more 

1,200 

9,800 

89 

1  or  more 

2,200 

11,000 

100 

the  number  of  threads  has  to  be  Increased  from  500  to  1,200.  Thus  after 
500  threads,  the  gain  in  covering  more  traffic  becomes  relatively  insig¬ 
nificant  even  when  a  large  number  of  threads  are  added.  It  is  therefore 
suggested  that  for  this  first  phase  of  study,  the  design  can  be  started 
with  pairs  exchanging  10  or  more  flights  (about  250  pairs),  then  go  on 
to  pairs  exchanging  5  or  more  flights,  and  so  on. 

Although  only  250-500  community  pairs  out  of  2,200  pairs  account 
for  a  major  portion  of  the  high  altitude  traffic,  these  are  still  quite 
large  numbers  to  be  considered  simultaneously.  A  possible  approach  to 
decompose  the  250-500  pairs  into  further  smaller  groups  is  discussed 
below. 

v  s 

A  study  of  the  requested  altitudes  by  flights  between  various  com¬ 
munities  has  indicated  that  those  community  pairs  which  are  less  than 
300-400j  miles  apart  seldom  use  altitudes  above  28,000  ft.  On  the  other 
hand,  community  pairs  more  than  about  400  miles  apart  seldom  request  f 
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altitudes  of  less  than  28,000  ft.  Therefore,  except  for  the  number  of 
route  intersections,  all  other  network  attributes  can  be  computed  inde¬ 
pendently  for  the  two  sets  of  community  pairs.  After  establishing  these 
attributes  separately  for  each  set,  the  two  sets  can  be  considered  si¬ 
multaneously  to  establish  only  the  total  number  of  route  intersections, 
and  then  the  other  attributes — already  computed  separately--can  be  com¬ 
bined  algebraically  to  obtain  a  total  configuration. 

A  preliminary  list  of  community  pairs  exchanging  10  or  more  flights 
per  day  (above  18,000  ft)  is  included  in  Appendix  H.  The  list  is  based 
on  peak  IFR  data  supplied  to  us  by  NAFEC.  Part  1  of  that  list  shows 
those  pairs  which  are  more  than  400  miles  apart  and  are  expected  to  use 
basically  altitudes  of  28,000  ft  or  above.  In  Part  2  are  listed  those 
pairs,  which  are  less  than  400  miles  apart  and  are  expected  to  utilize 
basically  the  altitudes  between  18,000  and  28,000  ft.  It  is  stressed 
that  the  abovementioned  groupings  were  made  considering  only  the  high 
altitude  flights  (18,000  ft  and.  above)  and  are  presented  as  a  possible 
preliminary  approach  to  decomposing  the  community  pairs  for  convenient 
implementation  of  the  methodology.  Further  study  of  the  traffic  patterns 
and  projected  traffic  demands  might  suggest  other  suitable  groupings. 

Assuming  that  the  suggested  approach  is  acceptable,  the  network  for 
the  community  pairs  exchanging  10  or  more  flights  can  be  established  in 
the  following  steps: 

Consider  the  125  community  pairs  of  Part  1  in  Appendix  H. 

Using  the  proposed  methodology  and  computer  program,  establish 
a  few  configurations  that  appear  to  have  the  most  desirable 
combinations  of  network  attributes.  Let  these  selected  con¬ 
figurations  be  called  A,  B,  C  ...  etc. 

Consider  the  community  pairs  in  Part  2.  To  those  pairs  in 
which  the  reference  community  was  also  a  reference  community 
in  Part  1,  assign  the  standard  headings  corresponding  to 
Configuration  A.  To  other  pairs,  assign  various  suitable 
reference  headings  and  establish  the  network  attributes.  Choose 
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that  combination  of  reference  headings  which  in  combination 

with  Configuration  A  gives  the  most  desirable  set  of  attributes 

Let  the  superposition  of  Configuration  A,  along  with,  the  routes 

* 

s  based  on  the  best  headings  from  set  2,  be  called  A  .  The  con- 
7  figuration  A  has  all  the  community  pairs  of  Part  1  and  2. 

(3)  Repeat  .2  except  that  this  time  assign  the  headings  of  Configura 
tion  B  to  the  reference  communities  common  to  both  sets. 

/  5|C 

.  Establish  the  configuration  B  . 

*  *  * 

(4)  After  establishing  the  configurations  A  ,  B  ,  C  ...  by  suc¬ 
cessively  applying  the  proposed  methodology  and  the  program, 
compare  these  configurations  and  select  the  one  which  appears 
to  have  the  ,best  combination  of  attributes  (and  possibly  has 
some  other  desirable  features  not  included  in  the  program). 


G.  Matching  a  Route  Network- and*  a  VORTAC  Grid 


The  route  network  design  methodology  considered  so  far  was  based  on 

the  assumption  that  adequate  VORTAC  coverage  is  available  everywhere  in 

the  United  States.  At  present,  there  exist  about  280  high  altitude 

.VORTACs  on  the  U.S.  mainland  and  indeed  these  provide  more  than  adequate 

coverage.  However,  as  discussed  in  Part  2  of  the  report,  adequate  cover- 
> 

age  above  18,000  ft  can  be  provided  with  as  few  as  100  VORTACs  if  these 
are  located  in  a  suitable,  orderly  pattern,  e.g., ,a  triangular  pattern. 


One  question  naturally  arises:  would  a  certain  orientation  of  a 
proposed  triangular  VORTAC  grid  support  a  given  route  network  better 
than  others?  To  answer  this  question,  it  is  necessary  to  establish  the 
route  supporting  characteristics  of  a  VORTAC  grid  for  various  orienta¬ 
tions  of  the  grid.  A  proposed  measure  of  the  route  supporting  character¬ 
istics  of  a  VORTAC  grid  for  a  given  route  network  is  the  percentage  of 
route  mileage  or  .flight  mileage  flyable  with  a  certain  accuracy,  e.g., 
with  a  certain  route  width.  Thus,  for  a  given  route  network,  if  a 
certain  orientation  of  a  proposed  triangular  VORTAC  grid  results  in 
higher  route  mileage  or  higher  flight  mileage  flyable  with  a  specified 
route  width,  e.g.,  8  nmi,  it  can  be  stated  that  this  orientation  is 
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relatively  superior  to  other  orientations  in  terms  of  route  supporting 
characteristics. 

In  Appendix  I  of  this  report,  simplified  relationships  have  been 
developed  by'  which  the  percentage  of  a  route  flyable  with  a  specified 
route  width  can  be  calculated  as  a  function  of  the  side  length  of  the 
triangular  VORTAC  grid  and  the  orientation  of  the  .route,  relative  to  a 
reference  row  of  the  triangular  grid. 

A  computer  subprogram,  based  on  the  relationships  shown  in  Appendix  I, 
has  been  developed  by  which  the  percentage  of  route  mileage  and  flight 
mileage  flyable  with  a  specified  route  width  is  computed  for  a  given 
•side  length  and  orientation  of  a  triangular  grid.  The  percentages  of 
route  and  flight  mileage  flyable  with  specified  accuracies  can  be  used 
as  a  measure  of  goodness  of  match  between  a  route  network  and  a  VORTAC 
grid.  •  ■ 

H .  Establishment  of  Way  Points  and  Changeover  Points 

After  deciding  upon  a  route  network  and  a  VORTAC  grid  on  the  basis 
of  broad  network  and  grid  attributes,  the  next  step  would  be  to  establish 
way  points  and  changeover  points  along  each  route.  The  development1 of  a 
detailed  computer  program  to  accomplish  this. task  is  beyond  the  scope  of 
this  project.  However,  the  NAFEC  team  has  already  developed  several 
mathematical  relationships  and  subroutines  necessary  for  this  task. 

Here  we  only  outline  the  steps  necessary  to  establish  the  way  points  and 
changeover  points  with  reference  to  NAFEC ' s  technical  notes  4  and  5 

j 

(listed  in  the  bibliography): 

(1)  Pick  a  route  or  route  section  i.  The  latitude  and  longitude 
of  its  start  and  end  points  are  known. 

(2)  Using  the  relationships^  of  NAFEC  technical  note  4,  find  the 
normal  distances  of  all  the  VORTACs  in  the  vicinity  of  the 

I 
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roufe  under  consideration.  The  latitude  and  longitude  of  all 
VORTACs  are  known  from  VORTAC  grid  program  output. 

(3)  Select  those  VORTACs  which  are  at  most  a  distance  K  nmi  away 
from  the  route,  so  as  to  pick  up  only  those  VORTACs  which  can 
provide  reasonable  coverage  to  the  route.  A  typical  value  of 
K  could  be  130  nmi. 

(4)  From  all  those  VORTACs  which  are  at  most  K  nmi  away,  pick  up 
the  best  set  (in  some  sense,  e.g.,  giving  minimum  number  of 
VORTACs).  The  projections  of  these  VORTACs  on  the  route  con¬ 
stitute  suitable  way  points  for  the  route. 

(5)  We  now  have  a  known  set  of  VORTACs  assured  to  support  the  route 
under  consideration  and  the  corresponding  way  points.  Using 
the  relationships  of  NAFEC  technical  note  5,  establish  the 
changeover  points  producing  equal  cross-course  error. 

The  above  outline  is  presented  as  a  preliminary  approach.  The  way  points 
and  the  changeover  points  so  established  will  have  to  be  flight  checked 
before  they  can  be  finalized. 
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VI  COMPUTER  PROGRAMS  TO  IMPLEMENT  THE  METHODOLOGY 

/ 

A.  Program  Section  1:  To  Establish  Intersection  Points  and  Angles 

Note:  for  the  derivation  of  the  various  formulas  and  relationships 

used  below,  refer  to  Appendix  F.  . 


1.  Input  Data 

•  The  geodetic  latitudes  and  longitudes  of  the  community  air¬ 
ports  are  input.  A  route  or  a  route  section  i  is  specified  by  four  numbers 


0 

0 

0 

e 


* 

il 

* 

11 
’* 

12 
* 

i2 


latitude. at  one  community 
longitude  at  the  same  community 
latitude  at  the  other  community 
longitude  at  the  other  community. 


i  =  1,  2,  3,  ...  n  routes. 


Convention:  always  choose  0  <  0  (i.e. .  9  always  to  the  east  o: 

il  i2  il 

* 

0  ).  Distance  D  in  nautical  miles  (terminal  area  radius). 

1  ^ 


Note:  The  program  converts  the  geodetic  latitudes  and  longitudes  to 

1  geocentric  system  for  calculating  intersection  points,  merging, 
demerging,  route  lengths,  and  so  on,  and  converts  them  back  to 
geodetic  units  before  outputting  the  results.  Thus  all  outputted 
latitudes  and  longitudes  are  geodetic. 


The  mathematical  relationships  described  here  are  based  on  east  to  west 
network  design  considerations.  However,  only  a  few  minor  changes  are 
necessary  to  make  these  relationships  applicable  to  west  to  east  net¬ 
work  design.  The  necessary  modifications  are  incorporated  in  the 
computer  program. 
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2. 


What  the  Program  is  Required  to  Do 


•  Calculate  the  latitude  and  longitude  (0  ,9  ) ,  (0  ,0  ) 

il  il  i2  i2 

of  the  start  and  end  points  of  each  direct  route;  a  distance  D  nmi  from 


the  community  airports  as  follows: 


First  calculate: 


*  *  *  *  *  * 
d  =  cos  [sin  0  sin  0  +  cos  0  cos  0  cos  (0  -  9  )1 

i  il  i2  il  i2  i2  il 


then  calculate: 


fil  = 


cos 


*  * 

, sin  0  -  sin  0  cos  d 

-11  12 _ _il _ i 

*  * 

cos  0  sin  d 
il  i 


f  i2  = 


,  -1 
cos 


*  * 

sin  0  -  sin  0  cos  d 

il  i2  i 


cos  0  sin  d 
i2  i 


then  assuming  1°  =*  60  nmi: 

0  =  sin  ~  |sin  0  cos(--) 

11  L  il  \60/ 

0  =  sin  1  sin  0  _  cos/77) 

12  L  i2  \  60  / 


=  sin  ^  sin  0  cos(--)  +  cos  0  sin/—  \  cos  Til] 

il  L  il  \60/  il  \60/  J 


*  /D  \  *  /D  \ 

sin  0  cos( — 1  4-  cos  0  sinl  —  J  cos  Ji 

\60  / 


i2 


*  -II 

0  =9  +  cos 

il  il 


cos 


(-) 

\60  / 


sin  0  sin  0 

il  il 


cos  0  cos  0 

il  il 


*  -1| 
9=0-  cos 
i2  i2 


cos 


(k)'~  sln  0i2  sin  0i: 


cos  0  cos  0 

i2  i2 


(1) 

(2) 


(3) 


(4) 
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•  Calculate: 


X  =  cot 

0 

cos 

0 

il 

il 

il 

y  =  cot 

0 

sin 

0 

il 

il 

il 

X  =  cot 

0 

cos 

e 

i2 

i2 

12 

y  =  cot 

•0 

sin 

9  • 

i2 

i2 

i2 

i  =  1,  2,  3,  . .  . ,  n 

For  i  =  1,  2,  3,  n  and  for  j  =  (i  +  1),  (i  +  2) , 

(i  +  3),  .  n,  establish  numbers  b 


b  =  1  if  Routes  i  and  j  intersect 
ij 

=  0  if  Routes  i  and  j  do  not  intersect 

For  Routes  i  and  j,  given  by: 

0  ,  9  ,  ^  ,  0 (Route  i) 

ll  ll  i2  i2 


0jl  ’  9jl  ’  0j2  ’  9j2  (R°Ute  J> 


•  Calculate: 

(x.Q  -  x  )(y  -  y  )  -  (x  -  x  )(y  -  y  ) 

j2  i2  j2  jl _ j2  jl  j2  i2 

ij  ”  (x  -  x  )(y  -  y  )  -  (x  -  x  )(y  -  y  ) 

ll  i2  j  2  jl  j2  jl  ll  i2 


(x  -  x  ) (y  -  y  )  -  (x  -  x  ) (y  -  y  ) 

il  12  j'2  yi2 _ j2  i2  il  i2 

ij  <x  *  x  o)(y-o  •  y.  >  -  <x  -  x  )(y.n  -  y,0) 

il  i2  j2  Jl  j2  jl  ll  i2 


M .  . 


Check  whether  0  <  X  <1 

ij 

and  0  <  (j.  <  1 

i  j 


If  so.  Routes  i  and  j  intersect;  otherwise,  they  do  not 
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Print  a  listing  of  intersecting  routes  in  a  suitable  format. 


For  all  Routes  i.  i  that  intersect,  find  the  latitude  0  and 

ij 

longitude  0  ,  as  follows: 

ij 


0 

ij 


0. 

IJ 


-1 

tan 


-1 

cot 


X  y 

+ 

<i 

-  X 

>y  o 

ij  il 

ij 

i2 

X  x 

+ 

(1 

-  X 

)  x 

ij  il 

ij 

i2 

X.  .y 

+ 

(1 

-  X 

)y .  _ 

ij  il 

ij 

i2 

sin  0 

ij 

Calculate  the  intersection  angle  &  of  intersecting  Routes 

ij 


i,  j  as  follows: 


First,  calculate: 


x  =  X  x  +  (1  -  X  )  x 
ij  ij  il  ij  ±2 


Y.  .  =  X.  .  y .  +  (1  -  X.  .)  y 

ij  iJ  il  ij  i2 


Then, 


Ot  -  cos 
IJ 


-1 


A  +  B  +  C 


/D  E 


Where 


A  - 

(y 

-  y.n>(y. 

-  y 

ij 

n  u 

ji 

B  = 

(x  . 

-  x  )  (x.  , 

“  x  ) 

ij 

ll  IJ 

Ji 

C  = 

(x  y 

.  -  y .  x.  . 

)  (x  y 

-  V  x  ) 

il 

ij  ii  ij 

Jl 

ij  jl 

ij 

D  = 

C(y  . 

-  y  )2  + 

(x.  .  - 

x  )2  + 

(x  y  .  ■ 

-  y.x.  ,)2] 

ij 

ii 

ij 

il 

ii  ij 

ii  ij 

E  = 

{(y.  . 

-  y  )2  + 

(X.  .  - 

X  )2  + 

(XV.  .  ■ 

■ 

< 

ro 

ij 

ji 

IJ 

jl 

jl  IJ 

Jl  iJ 
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Print  a  list  of  intersection  coordinates  in  a  suitable  format. 
Print  a  list  of  intersection  angles  in  a  suitable  format. 

B.  Program  Section  2:  To  Establish  Network  Attributes 

Note:  for  the  derivation  of  the  various  formulas  and  expressions 

used  below ,  refer  to  Appendices  C  and  D. 

1.  Input  Data 

Matrix  A  (n  x  5) 

a  =  number  of  flights  on  Route  i  at  altitude  L  during  h  hours 
iL 

) 

L  =  1,  2j  3,  4j  5. 


Listing  of  intersection  routes  and  intersection  angles  from 
Program  Section  1. 

Constants: 

Constant  _ Description _ _  Typical  Value 


r 

V 

h 

k 


a 


b 


Radar  separation  minima  in  miles  5 

Average  speed  of  aircraft  in  miles  per  hour  500 

Duration  of  traffic  count  in  hours  5 

Cost  per  potential  conflict  in  dollars  per 

conflict  10 

Allowable  conflict  duration  in  fraction  of  an 

hour  0.01 

Allowable  number  of  conflicts  during  h  hours  5 
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2. 


What  the  Program  is  Required  to  Do 


The  program  calculates  and  prints  the  information  listed  below. 


a.  Number  of  Potential  Conflicts  at  Various  Intersections 


•  The  matrices  E  (n  x  n) (L  .=  1,2, 3, 4,5)  with 
L 


a 


2r  /  ij  \ 

c  ,  L  =  —  b  ( sec  -  a  a 

ijJ  hV  ij  2  iL  jL 


•  The  sum  of  the  rows  of  E  (E  e  .L  is  the  total  number 

L  j  ij 

of  conflicts  along  route  i  at  Level  L). 

•  The  sum  of  the  conflicts  e  =  1/2  Y.  £  e  ,L  (the  con- 

L  i  J  ij 

flicts  at  Level  L) . 


b.  Average  Conflict  Duration 


O' 


r  i  j 

•  The  matrix  F  (n  X  n)  with  f  =  15  tt  —  b  cosec  - 

ij  V  ij  2 


c.  Conflict  Duration  Time 


•  The  n  X  n  matrices  G  (L  =  1,2, 3, 4, 5)  with 

L 


60  tt  r 

g  ,  L  =  - - —  b  (cosec  o  )a  a 

ij  .  .2  ij  ij  iL  JL 

hV 


The  sum  of  the  rows  of  G  (£  g  ,L  is  the  total  con- 

L  j  ij 

flict  duration  along  Route  i) . 


The  sum  of  total  conflict  time  at  Level  L? 


g.  =  1/2  ?  ?  g.  ,}L 
L  1  J  1 J 
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d. 


Intersection  Capacity  Based  on  Conflict  Duration,' 

•  The  matrix  H  (n  x  n). 

2  '  ,‘ 

a  V 

h.  —  —  —  b  sin  O' 
ij  TT  2  ij  ij 

e.  Intersection  Capacity  Based  on  Conflict  Numbers 


•  The  matrix  M  (n  x  n), 

a 

b  v  .  ij 

m  =  —  —  b  cos  - 

ij  2h  r  ij  2 

f .  Total  Flight  Mileage 

•  The  value  for  E  f  j l 


ill 


where 


f  =  £  a 
i  L  iL 


and 


J l  =  3440  cos  ^  [sin  0  sin  0 
i  il  i2 


+  COS  0.  COS  0'  cos  (0  -  8  ) 

il  i2  i2  il 


C.  Program  Section  3:  Route  Merging 

Note:  for  the  derivation  of  the  various  formulas  and  expressions 

used  below5  refer  to  Appendix  B. 


1.  Input  Data 

The  input  data  consist  of  the  following: 


•  Input  data  as  in  Program  Section  1  and  specification  of  a 
reference  heading  or  a  reference  route  for  each  community. 

•  Constant  8;  the  departing  angle  of  demerging  routes. 

(«  10°) 


A 

heading  increment^  degrees  (note:  make  sure  that  9  ^  — ) 

2 
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2.  What  the  Program  is  Required  to  Do  , 

•  Let  the  reference  direction  of  a  community  be  specified  as 
Route  j.  . 

,  Calculate  the  heading  y 

j 

(the  angle  between  north  and  the  route;  measured  in  the  clockwise 
direction): 


First,  calculate  the  length  d  . 


d  =  cos  Tsin  0  sin  0  -h  cos  0  cos  0  cos  (9  -  9  )] 

j  jl  J2  .  Jl  J2  j2  jl 


Then, 


-1 

cos 


sin  0  -  sin  0  cos  d 

_ jl _ j 

cos  0  sin  d 

jl  j 


If  the  reference  direction  of  a  community  is  already  specified 
as  a  heading  j,  the  above  calculations  are  not  necessary. 


•  Define  standard  heading  directions  h  (n  =  0, 1 , 2, 3, .  .  .  , ) : 

j.n 


h  =  (y  +  A  n)  if  Y  +  An<  360° 
J>n  J  j 


=  (y  +  A  n  -  360)  ^  if  y  +  A  n  >  360° 

j  j 


All  other  routes  from  this  community  will  be  partially  merged 
along  one  or  the  other  heading  directions,  as  shown  in  Figure  3  and 
described  below. 

,  / 

A  Route  k  will  be  merged  with  the  heading  direction  for  which 

h_  -  y  <  0  .  Let  the  latitude  and  longitude  of  the  point  from 

j,n  k 

which  route  k  separates  be  designated  as  0  ,0  .We  want  0 

jn,k  -jn,  k  *  jn,k 

and  0  'in  terms  of  0  ,0  and  v  and  v 

.  jn, k  jl  >  jl  j  ^k 
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7.  +  20  n 

(0.  6.  .)  ' 

^jn,k  jn,k 


SA-1096-11 

FIGURE  3  DEMERGING  POINT  OF  A  MERGED  ROUTE 


Firsts  calculate: 

g  =  d  [cos  (v  -  h  ,  )  -  [ | sin  (V  -  h.  )|]  cot  0} 

k  k  k  j^n  _l  k  n  1 

-■“T  [ I  sin  (y  -  h .  ) | }  cot  9 

60  k  j;n  1 

where 

-  d  =  cos  1  [sin  0  sin  0  +  cos  0  cos  0  cos  (0,  -  9  )1 

k  L  kl  k2  kl  k2  k2  kl  J 

(already  calculated).  Note  that  g  is  the  length  of  the  merged  part. 

k 
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If  g  >  0,  then 
k 

0  =  sin  [ sin  0  cos  g  +  cos  0  sin  g  cos  (360  -  h  )] 

jn,k  jl  k  jl  k  j,n 


0 

3n> 


k 


0  +  cos 

jl 


cos  g  -  sin  0  sin  0 

_ k _ jn,k _ jl 

cos  0  cos  0 

jn,k  jl 


If 


0, 


then  the  corresponding  route  is  left  unmerged. 


nated  as 


Let  all  routes  that  are  mergeable  with  heading 
k,  m  . . . ,  such  that 


be 


redesig- 


e  •  <  e  <  e 

jn,k  jn,£  jn,m 


Then^  the  original  routes,  k,  SL i  m,  from  the  conununity  under  consideration 
will  be  redirected  as  shown  in  Figure  4. 


SA- 1096-5 


FIGURE  4  REDIRECTIONS  OF  ORIGINAL  ROUTES 
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Route  Section  Between 


With  Flows 


2  a 

sL 


(VV 

(0jnJk^jn3k) 

S  =  £, 

»i_i 

m,n,  o 

^jnjk^jn.k5 

(0k2’6k2) 

a 

kL 

(0jn,k"0jnJk> 

<0jn>A,Sjn,l) 

(2a  ) 

sL 

-  a 

kL 

(0jn,i'9  jn-J£) 

(W 

*1 L 

(0\  ,,e.  .) 

jn,£  jn  ,1 

(0.  ,6.  ) 
jn,m  jn, m 

(2a  ) 

sL 

-  a 

kL 

(0.  ,e.  ) 

jn,m  jn,ni 

«„2’V 

a 

mL 

(0.  ,e .  ) 

jn,m  jn}m 

‘VV' 

a 

jL 

In  the  original  table  of  routes,  the  routes  emanating  from  the 
community  under  consideration  are  modified  as  above.  The  new  intersecting 
points  are  calculated  by  using  Program  Section  1.  The,  new  network  at¬ 
tributes  are  calculated  by  using  Program  Section  2.  The  results  are 
tabulated  for  study.  * 

The  entire  process  of  the  route  merging  program  is  continued 
by  successively  selecting  different  communities  and  assigning  suitable 
reference  directions  to  each  community,  establishing  the  corresponding 
standard  headings,  and  then  merging  the  routes  from  the  community  ac¬ 
cording  to  the  above  method  until  all  communities  are  taken  into  account. 

The  designer  can  select  various  combinations  of  reference  routes 
or  headings  and  study  the  resulting  network,  compare  various  network 
attributes,  and  then  select  those  configurations  which  appear  to  have  the 
most  desirable  combination  of  attributes. 
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VII 


CONCLUDING  REMARKS 


The  methodology  and  the  computer  program  developed  so  far  are  of  a 
preliminary  nature  and  are  based  on  some  simplifying  assumptions.  The 
purpose  of  this  section  is  to  review  briefly  some  of  the  significant 
simplifying  assumptions  and  discuss  their  effects  on  the  validity  of  the 
design  methodology. 

Following  is  a  list  of  some  important  simplifying  assumptions. 

•  The  aircraft  flows  and  speeds  have  been  represented  as  average 
flows  and  speeds. 

•  The  effects  of  controller  judgmental  factors  have  not  been 
taken  into  account. 

•  Restricted  areas  have  not  been  considered  explicitly. 

•  Eastward  and  westward  route  networks  are  considered  independently. 

•  The  merged  routes  generated  by  the  present  methodology  occasion¬ 
ally  intersect  with  each  other  twice.  Such  situations  have  to 
be  rectified  manually. 

We  now  discuss  each  of  the  above  assumptions  in  the  order  stated. 

A .  Average  Flows  and  Speeds 

If  a  network  were  to  be  developed  based  on  air  traffic  as  it  is 
actually  distributed  in  time  and  space,  one  would  either  have  to  employ 
elaborate  simulation  techniques  or  develop  statistical  models  requiring 
extensive  traffic  data.  Both  these  approaches  were  considered  beyond 
the  scope  of  this  project.  A  proposed  approach  to  account  for  dynamic 
variations  in  traffic  flow  is  to  use  "peak  hour”  flows  as  explained  in 
Appendix  G,  producing  a  "conservative”  network  based  on  the  assumptions 

that  peak  flows  occur  during  the  same  hour.  However,  since  in  reality 

f 


Preceding  page  blank 
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this  will  not  be  the  case,  it  will  be  advisable  to  reexamine  the  designed 
"conservative"  network  to  ensure  that  no  unwarranted  steps  were  taken. 

It  is  suggested  that  sample  routes  designed  using  the  proposed  methodology 
should  be  tested  in  the  NAFEC  simulation  facilities.  Such  testing  will 
provide  valuable  insight  and  will  be  helpful  in  establishing  the  areas 
of  potential  improvement. 

B .  Controller  Judgment 

•  ATC  controllers  play  a  significant  role  in  shaping  the  traffic  flow 

patterns.  However,  some  of  the  actions  of  ATC  controllers — particularly 

the  ones  based  on  their  judgment  and  experience — are  extremely  difficult, 

if  not  impossible,  to  model  in  terms  of  analytical  expressions.  As  such, 

human  factors  were  not  included  in  the  preliminary  design  methodology. 

* 

In  a  related  project,  SRI  is  currently  developing  techniques  to  estab¬ 
lish  the  constraints  on  the  ATC  capacity  due  to  human  controllers.  It 
is  conceivable  that  some  of  these  techniques  could  be  incorporated  in 
the  route  network  design  methodology  in  a  separate  study. 

C.  Restricted  Areas 

In  the  early  stages  of  this  project,  it  was  intended  that  restricted 
areas  would  be  taken  into  account  explicitly.  However,  it  now  appears 
that  the  logic  and  the  computer  storage  requirements  to  accomplish  the 
bending  of  routes  around  restricted  areas  will  be  too  extensive  to 
justify  the  cost.  It  is  suggested  that  for  the  preliminary  design,  the 
bending  of  routes  around  the  restricted  areas  be  implemented  manually, 
i.e.,  first  design  the  routes  without  considering  the  restricted  areas 


SRI  Project  8181,  A  Methodology  for  Evaluating  the  Capacity  of  ATC 
Systems.  * 
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and  then  manually  modify  those  routes  which  are  crossing  restricted  areas. 
For  example,  suitable  route  sections  around  the  restricted  area  would  be7 
selected  manually  and  then  the  network' attributes  would  be  tested,  using 
the  program  to  select  a  proper  modification. 

D.  Route  Directions 

The  present  program  is  meant  to  design  the  eastward  and  westward 
networks  independently.  The  rationale  for  this  approach  has  already  been 
discussed  in  Section  V-E. 

Here  we  wish  to  point  out  that  the  two  versions  of  the  program  can 
be  used  serially  to  merge  half  lengths  of  the  routes  using  the  eastward 
version  and  the  remaining  half  lengths  using  the  westward  version.  For 
example,  the  direct  route  between  community  pairs  New  York  and  San  Francisco 
could  be  divided  into  two  halves.  The  part  from  New  York  to  the  middle 
point  can  be  merged  partially  along  a  standard  heading  from  New  York 
(westward  direction),  and  the  part  from  San  Francisco  to  the  middle  point 
can  be  merged  partially  along  a  standard  heading  from  San  Francisco 
(eastward  direction).  Note  that  the  complete  route  between  New  York  and 
San  Francisco  in  this  approach  will  consist  of  four  sections,  i.e., 
three  zigs.  It  will  be  useful  to  apply  this  approach  to  a  sample  network 
and  to  study  the  resulting  network  for  comparison  purposes. 

E .  Merged  Route  Intersections 

In  the  present  methodology  and  the  associated  computer  program, 
routes  are  partially  merged  along  standard  headings,  thereby  creating  a 
bend  in  the  route.  Thus,  a  route  that  originally  consisted  of  a  single 
straight  line  is  changed  into  a  route  with  two  straight  line  sections. 
Occasionally,  two  routes  i  and  j  that  originally  crossed  each  other  only 
once,  or  did  not  cross  at  all,  may  cross  each  other  twice  because  of  the 
bends  in  either  or  both  of  them,  as  illustrated  in  Figure  5. 
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FIGURE  5  CROSSING  OF  A  ROUTE  TWICE  WITH  ANOTHER  ROUTE  BECAUSE 
OF  THE  BEND 


In  this  figure,  routes  i  and  j  were  not  crossing  originally  when 
direct  path  was  provided.  When  route  i  is  merged  partially  along  some 
standard  heading  resulting  in  the  route  shown  with  broken  lines,  two 
intersections,  A  and  B,  result.  This  obviously  is  an  undesirable^  situ¬ 
ation.  Rectification  of  such  situations  through  computer  programs  would 

\ 

require  a  rather  complicated  logic  and  does  not  appear  to  be  cost  bene¬ 
ficial  at  this  stage.  It  is  suggested  that  for  preliminary  design  pur¬ 
poses,  these  situations  be  manually  corrected  by  merging  the  part  between 
A  and  B  along  the  route  j.  To  assist  the  designer  in  locating  such  un¬ 
desirable  situations,  a  provision  has  been  made  in  the  program  to  indi¬ 
cate  those  sets  of  intersections  and  routes  which  must  be  rectified.  A 
further  check  can  be  made  by  studying  the  plot  of  the  network.  After 
studying  the  list  of  all  routes  and  intersections  requiring  modifications, 
the  designer  may  manually  specify  appropriate  modifications. 

Although  the  developed  methodology  is  not  comprehensive  and  is  based 
on  several  simplifying  assumptions,  it  is  the  first  systematic  computer 
based  methodology  for  designing  large  route  networks.  The  methodology 

is  modular  in  structure  so  that  various  parts  can  be  changed  or  improved  t 

independently,  e.g.,  if  an  improved  conflict  model  is  later  developed, 

it  can  easily  replace  the  original  one  without  requiring  any  major  change 

i 

in  the  program. 
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Appendix  A 


PROBLEM  DEFINITION  AND  GROUND  RULES 

This  appendix  describes  a  tentative  set  of  ground  rules  for  the 
development  of  the  methodologies  for  Task  3.  The  rules  were  developed 
during  the  week  beginning  26  April  1971  and  reflect  agreement  among 
the  NAFEC  and  SRI  project  team  members  on  problem  definition — the 
requisite  starting  point  for  the  development  of  a  systematic  and  logical 
network  design  methodology. 

1 .  Traffic  Data  Base 

Resolution :  The  IFR  peak  day  traffic  tapes  will  be  used  as  the 

data  base  for  representing  traffic  demands. 

Discussion :  The  IFR  peak  day  traffic  tapes  were  chosen  to  provide 

a  data  base  because  they : 

•  Are  frequently  used  for  characterizing  traffic  levels  at 
Air  Route  Traffic  Control  Centers  ( ARTCC )  and  Airport 
Traffic  Control  Towers  (ATCT) . 

•  Will  provide  a  conservative  (or  worst-case)  system  design. 

•  Include  traffic  from  all  three  classes:  air  carrier, 
military  aviation,  and  general  aviation. 

•  Contain  most  of  the  required  traffic  demand  parameters  at 
the  community  centers  and  are  in  a  form  that  is  easily 
adaptable  for  the  present  purposes. 

•  Can  readily  be  scaled  to  represent  forecasted  traffic  demands. 

The  community  centers  and  community  pairs  to  be  included  in  the 
design  of  the  route  network  should  be  determined  on  the  basis  of  traffic 
level.  This  traffic  demand  specification  has  been  incorporated  into 
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the  traffic  demand  model  by  rank  ordering  the  community  centers  and 
community  pairs  according  to  IFR  high-altitude  (18,000  ft  and  above) 
traffic  levels. 

2 .  Interface  with  Terminal  Areas 

a .  Departure /Arrival  Points 

Resolution :  Departure/arrival  points  will  be  located  at  a 

radius  of  50  nmi  from  the  geometrical  center  of  a  community.  Latitudes 
and  longitudes  for  a  set  of  discrete  departure/arrival  points  will  be 
derived  and  furnished  for  each  community  center  in  the  data  base. 

* 

Discussion :  NAFEC  will  develop  these  points  on  the  basis 

* 

of  traffic  density  and  the  angles  to  the  interfacing  communities. 
Departure/arrival  points  derived  in  this  manner  will  reflect  a  best 
choice  from  the  aircraft  user's  viewpoint.  In  practice,  geographical 
or  wind  factors  may  not  permit  this  choice.  Such  factors  will  be  con¬ 
sidered  in  later  phases  of  this  study. 

b .  Climb/pescent  Areas 

Resolution :  The  problem  of  climb  and  descent  is  to  be  con¬ 

sidered  primarily  as  a  terminal  problem;  it  will  not  be  handled  ex¬ 
plicitly  in  this  phase  of  the  design  of  the  route  network. 

Discussion :  Since  it  is  known  that  climb  to  and  descent  from 
the  high-altitude  airspace  occur  beyond  the  50-nmi  radius  that  is  being 
assumed  for  the  terminal  boundary,  provisions  will  be  made  for  this  inter¬ 
face.  One  possible  means  of  accomplishing  this  is  to  allot  a  certain 


The  rule  for  choosing  a  terminal’s  departure/arrival  points  (on  a 
circle  of  radius  50  nmi),  which  is  based  upon  predominant  (or  prime) 
traffic  flows,  is  described  in  NAFEC's  Technical  Note  No.  8. 
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amount  of  airspace  adjacent  to  the  routes,  for  an  additional  number  of 
miles  from  each  departure/arrival  point.  This  provision  is  a  reasonable 
approximation  for  the  purposes  of  this  phase  of  the  study,  and  it  should 
provide  sufficient  flexibility  when  the  problem  of  climb  and  descent  in 
the  terminal  areas  is  considered  explicitly  in  subsequent  phases  of  the 
design . 

3 .  Restricted  Airspace 

Resolution :  Restricted  airspace  will  be  considered  in  the  design 
of  the  route  network. 

Discussion :  Since  restricted  airspace  imposes  spatial  constraints 
(latitude,  longitude,  and  elevation)  on  the  location  of  routes,  it  is 
an  important  element  that  must  be  incorporated  in  the  design  of  the 
route  network  structure. 

The  data  that  define  restricted  regions  in  terms  of  their  area 
projection  on  the  ground  and  their  altitude  band  are  readily  available 
from  navigational  charts  and  other  sources.  These  data  are  now  being 
compiled  for  regions  in  the  high-altitude  (18,000  ft  and  above)  air¬ 
space.  For  the  purposes  of  high-altitude  route  network  design,  the 
ground  projection  for  each  restricted  region  can  be  represented  by  a 
convex  boundary.  This  convex  boundary  is  obtained  by  connecting  a  sub¬ 
set  of  ground  projection’s  vertices  with  great  circles,  as  illustrated 
in  Figure  A-l.  It  should  be  noted  that  this  entails  no  loss  of  gener¬ 
ality  for  the  case  of  high-altitude  routes,  given  the  geometry  and  size 
of  actual  restricted  regions. 

4.  Route  Widths 

Resolution1:  Route  widths  will  be  assumed  to  be  constant  for  the 
design  of  the  route  network  . 
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FIGURE  A-1  GROUND  PROJECTION  FOR  A  RESTRICTED  REGION 

Discussion :  In  practice,  route  width  is  determined  by  the  geo¬ 

metrical  relationship  of  the  supporting  VORTACs  to  the  particular  route 
Since  the  design  of  the  route  network  in  Task  3  is  to  be  independent  of 
VORTAC  geometry,  it  will  be  assumed  that  the  route  widths  are  some  con¬ 
stant  value  representative  of  typical  ATC  separation  requirements  (e.g. 
8  to  12  nmi).  Furthermore,  in  the  Task  3  methodology,  route  widths  can 
be  varied  parametrically  in  the  design  algorithm-  The  relationship  of 
route  widths  to  VORTAC  geometry  will  be  considered  explicitly  in  the 
Task  4  methodology. 

5 .  Altitude  Utilization  Rules 

Resolution :  The  hemisphere  rule  will  be  used  as  a  starting  point 

from  which  other  altitude  utilization  rules  will  be  investigated  for 
possible  inclusion  in  the  Task  3  methodology. 

Discussion :  In  the  development  of  the  Task  3  methodology,  it  is 

planned  to  investigate  other  altitude  utilization  rules  that  are  poten¬ 
tially  more  efficient  than  the  hemisphere  rule  in  the  use  of  available 
airspace.  That  is,  by  developing  an  algorithm  that  more  efficiently 
utilizes  the  available  altitudes,  one  could  expect  a  reduction  in  the 


54 


number  of  intersections  between  routes  or  in  the  expected  number  of  po¬ 
tential  conflicts  between  aircraft.  Such  a  reduction,  in  turn,  would 
lead  to  increased  capacity  for  the  resulting  route  network.  In  this 
approach,  the  hemisphere  rule  would  serve  as  a  basis  for  comparing  other 
altitude  utilization  rules. 

6 .  Atmospheric  Conditions 

Resolution :  Atmospheric  conditions,  i.e.,  weather,  wind,  and  jet 
stream,  will  not  be  considered  explicitly  in  the  initial  design  of  the 
route  network. 

Discussion :  Too  little  definitive  information  is  available  to  allow 
design  to  be  based  explicitly  on  atmospheric  conditions.  It  is  felt 
that  the  route  network  structure  to  be  developed  will  be  sufficiently 
flexible  to  enable  the  system  to  cope  with  atmospheric  conditions,  just 
as  the  present  system  does.  An  examination  will  be  made  in  later  phases 
of  project  work  to  determine  whether  this  capability  is  implicit  in  the 
design . 

7 .  Form  of  RNAV 

Resolution :  All  aircraft  will  be  considered  to  be  equipped  with 
rho-theta  type  RNAV  equipment  with  slant-range  correction. 

Discussion :  The  most  accurate  present-day  RNAV  is  the  rho-rho 

system.  The  implementing  of  PVORTAC  will  enable  rho-theta  type  RNAV 
to  match  the  accuracy  of  the  rho-rho  type  RNAV.  Further,  the  rho-rho 
system  requires  on  the  order  of  twice  the  ground  facilities  needed  in 
a  rho-theta  system.  Therefore,  the  trend  is  toward  the  latter,  and  it 
is  felt  that  this  is  the  RNAV  of  the  future.  Slant-range  correction 
will  certainly  be  in  all  future  systems. 
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8 .  VORTAC  Coverage 


Resolution :  It  will  be  assumed  that  full  VORTAC  coverage  is 

available  throughout  the  United  States. 

Discussion :  The  total  number  of  VORTACs  presently  available  in  the 

United  States  is  of  the  order  of  1100  of  which  300  VORTACs  are  meant 
for  high  altitude  routes.  With  their  present  locations,  the  existing 
VORTACs  provide  more  than  adequate  coverage  for  all  practical  purposes. 
The  justification  of  the  above  assumption  is  thus  obvious. 

9 .  High-Alt i tude  Level 

Resolution :  Only  air  routes  altitude  above  18,000  ft  will  be 
considered . 

Discussion :  It  is  planned  by  FAA  that  RNAV  routes  will  first  be 

introduced  for  high  altitude  routes.  The  altitude  of  18,000  ft  was 
chosen  as  a  preliminary  ceiling  to  define  a  suitable  finite  traffic  data 
base . 
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Appendix  B 

GROUPING  OF  THREADS  ON  THE  BASIS  OF  A  MINIMUM  DEPARTING 
ANGLE  BETWEEN  MERGED  AND  PRIMARY  THREADS 

The  purpose  of  this  appendix  is  to  develop  some  useful  relationships 
for  the  case  when  a  minimum  departing  angle  between  a  merged  and  a 
primary  thread  is  specified.  It  is  shown  that,  for  a  minimum  departing 
angle  of  10° ,  the  increase  in  the  length  of  the  thread  being  merged 
does  not  exceed  0.4  percent  of  the  direct  shortest  length  and  that,  for 
a  departing  angle  of  15°,  the  additional  length  is  less  than  0.9  percent. 

1 .  Some  Basic  Relationships 

The  basic  relationships  in  the  merging  of  threads  are  first  es¬ 
tablished  by  using  plane  trigonometry.  Later,  the  results  are  converted 
to  a  geocentric  system.  Referring  to  Figure  B-l,  let 

OP  be  a  primary  thread 

OA  a  candidate  thread  of  length  & 

Angle  AOP  =  O' 

Angle  G  =  minimum  specified  departing  angle 


FIGURE  B-1  A  PRIMARY  THREAD  (OP)  AND  A  CANDIDATE  THREAD  FOR  MERGING  (OA) 


r  Preceding  page&lankj 
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OC  =  ,  length  merged  along  the  primary  thread 

CA  =  ,  length  between  the  takeoff  point  and  Destination  A. 

AB  is  perpendicular  from  A  on  OP,  of  length  y  . 

Length  CB  =  x  . 

With  the  above  definitions,  the  following  relationships  can  immediately 
be  established: 


y  sin  9  =  y  =  £  sin  Of 
yl  2 


sin__Q- 

V  =  f  7 

^1  sin  8 


tan  9 


^2  &  sin  O' 


X  x  cos  Of  -  X 

£  1 


.*«  x  =  &  cos  O'  -  i  sin  Of  cot  9 


cos  of  -  sin  Qf  cot  9 


(B-l) 


(B-2) 


(B-3) 


(B-4) 


Equation  (B~4)  gives  the  length  of  the  part  of  the  thread  that  is  merged 
along  the  primary  thread  as  a  fraction  of  the  original  length. 


Let  us  now  calculate  the  increase  in  the  length  of  the  thread  when 
it  is  partially  merged.  We  have 


h  =  increase  in  the  length  =  +•  -  i 


sin"  o 


i-  cos  Of  -  &  sin  of  cot  9  +  Z  — - - r  -  & 


sin  9 


A  =  l 


(• 


A  sin  Of 
cos  Of  -  sin  Of  cot  9  +  — - - r  “  1 


sin  9 


(B-5) 
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Equation  (B-5)  gives  the  increase  in  the  length  of  Thread  OA. 

Dividing  (B-5)  by  l,  we  have 

cos  o'  -  sin  ex  cot  9  -  - -  -  1  )  .  (B-6) 

sin  9  I 

The  right-hand  side  of  (B-6)  can  be  simplified  by  using  standard  trigono¬ 
metric  relations: 


(B-7  ) 


Equation  (B-7)  gives  the  increment  in  the  length  of  the  thread  as  a 
fraction  of  the  original  length. 

X1  A 

In  Table  B~l,  values  of  “  X  100  and  “  X  100  have  been  calculated. 

-v  X> 

for  three  fixed  values  of  9:5° ,  10°,  and  15°. 


Table  B-l 


RELATION  BETWEEN  O',  x±/ ' H ,  AND  L/& 

FOR  FIXED  VALUES  OF  9,  THE  DEPARTING  ANGLE 


9  =  5° 

0  = 

10° 

e 

=  15° 

O' 

A 

i°o- 

X 

100r 

A 

1°0- 

X 

100r 

100A/4 

lOOx  A 

r 

0.5 

0.034 

90 

0.073 

95 

0.11 

97 

1 

.061 

80 

.14 

90 

.21 

93 

2 

.092 

60 

.24 

80 

.40 

87 

3 

.091 

40 

.32 

.55 

80 

5 

0 

0 

.38 

.77 

67 

7 

.32 

.86  . 

54 

9 

.14 

.83 

40 

12 

0 

.55 

20 

15 

- 

0 

0 
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In  Figure  B-2,  the  relationships  for  0  =  1CT  and  15°  are  shown 
graphically  for  further  clarity. 

From  Table  B-l  and  Figure  B-2  it  is  clear  that,  if  a  departing 
angle  as  large  as  15°  is  specified,  the  maximum  increase  in  the  length 
of  the  thread  being  merged  will  not  exceed  0.9  percent  (when  the  thread 
makes  an  angle  of  about  7°  with  the  primary  thread)  and  that,  if  the 
acceptable  departing  angle  is  taken  as  10° ,  then  the  increase  in  length 
will  not  be  more  than  0.4  percent  (when  the  thread  angle  is  5°). 


FIGURE  B-2  RELATIONSHIPS  BETWEEN  a,  A/Z,  AND  x,/« 
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2 .  The  Merged  Length  and  the  Coordinates  of  the  Take-off  Point  in 

Spherical  Coordinates 

Without  any  significant  loss  of  accuracy,  it  can  be  assumed  that 
Equations  (B“4)  and  (B~7)  are  valid  also  when  Points  0,  A,  P,  and  C  lie 
within  a  finite  area  of  the  earth's  surface.  Let  the  latitudes  and  longi¬ 
tudes  of  0,  A,  P,  and  C  be  (eA,0n),  (9  ,0  ),  (0o>0o)>  and  (9.,0.)» 
respectively,  as  shown  in  Figure  B-3 .  Let  N  be  the  north  pole.  Let 

V  and  y  be  the  headings  (the  angles  with  reference  to  local  north) 

1  2 

of  Routes  OA  and  OP.  Let  and  d^  be  the  great-circle  arc  angles  OA 
and  OP.  Then,  using  standard  relationships  of  spherical  trigonometry, 
we  have  (see  Figure  B-3)  : 


cos  3 [sin  0^  sin  0^  +  cos  0^  cos  0^  cos  (9^  -  6^)]  degrees  (B-8) 


d^  =  cos  ^Tsin  0  sin  0  +  cos  0  cos  0_  cos  (0 

2  L  0  2  0  2  2 


6^)]  degrees  (B-9) 


and  then 


Y1  =  cos 


-1 


sin  0^  -  sin  0Q  cos  d^ 


cos  0_  sin  d 
0  2 


degrees 


(B-10) 


Y  =  cos 


sin  0„  -  sin  0  cos  d 
2  0  2 


cos  0  sin  d 
0  2 


degrees 


(B-ll) 


Now,  Angle  or  =  <  AOP  =  (y^  “  y^) 


Substituting  d^  for  l  and  (y^  -  y^)  for  &  in  Eq .  (B-4),  we  get 


0C  =  x  =  d  [cos  <Y  “  Y  )  "  sin  (y  -  y  )  cot  ©]  degrees.  (B-12) 

^  -L  1  ^  1  Z 
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FIGURE  B-3  MERGED  LENGTH  IN  SPHERICAL  COORDINATES 

Considering  now  the  spherical  triangle  OCN,  we  have: 


cos  (90  -  0  )  =  cos  (x  )  cos  (90  -  0  )  +  sin  x  sin  (90  -  0  )  cos  Y0  (B-13) 

XX  U  X  U  £ 


sin  0.  =  cos  x  sin  0  +  sin  x  cos  0  cos  y 

i  10  10  2 


0  =  sin 

i 


-1 


sin  0  cos  x  +  sin  x  cos  0^  cos  y„ 
0  1  10  2 


degrees.  (B-14) 


Also 


cos  x  =  cos  (90  -  0  )  cos  (90  -  0  ) 
1.0  i 


+  sin  (90  -  0  )  sin  (90  -  0.)  cos  (0  -  0rt)  .  (B-15) 

O  l  i  0 
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Rearranging,  we  get 


9 

i 


=  0  +  cos 

0 


„  .  cos  x  -  sin  sin  0 

-l  I  l _ o _  i 

cos  0  cos  0 
0  i 


(B-16) 


Summarizing  the  results:  Equation  (B-12)  gives  the  length  of  the  merged 
part  of  a  thread  and  Equations  (B-14)  and  (B-16)  give  the  latitude  and 
longitude  of  the  take-off  point. 
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Appendix  C 

CONFLICT  PREDICTION  MODELS 

Under  radar  surveillance  conditions,  ATC  is  required  to  ensure  that 
two  aircraft  are  separated  at  least  by  5  nrai .  Whenever  this  condition 
is  violated,  the  two  aircraft  are  considered  to  be  in  conflict.  By 
using  this  definition,  simple  mathematical  models  for  predicting  the 
expected  number  of  conflicts  and  the  conflict  durations  at  an  intersec¬ 
tion  can  be  developed,  as  shown  below. 

1 .  Development  of  the  Mathematical  Model  for  Expected  Number  of  Poten¬ 
tial  Conflicts 

Let  : 

Ot  =  angle  between  Routes  AB  and  CD  (see  Figure  C-l). 

V  =  average  velocity  of  aircraft  on  the  routes  (nmi  per 
hour) . 

f^  =  expected  flow  rate  along  Route  AB  (aircraft  per  hour). 

f  =  expected  flow  rate  along  Route  CD  (aircraft  per  hour). 

X  =  minimum  separation  allowed  between  two  aircraft 
(typically,  5  nmi). 

Referring  to  Figure  C-l,  let  Aircraft  b  be  at  Intersection  0  and  flying 
toward  D.  We  ask  ourselves  where  should  Aircraft  a  not  be  along  Route 
AB  at  this  moment  so  that  (1)  a  conflict  between  Aircraft  a  and  b  is 
not  occurring,  (2)  a  conflict  will  not  develop  as  Aircraft  a  and  b 
continue  to  move  with  velocity  V  along  AB  and  CD,  respectively,  and 


FIGURE  C-1  TWO  INTERSECTING  ROUTES  FOR  EVALUATION  OF  POTENTIAL  CONFLICTS 

(3)  a  conflict  would  not  have  occurred  in  the  past  if  Aircraft  b  and  a 
were  moved  backwards.  We  denote  such  portions  as  critical  length  along 
AB. 


■  Consider  first  the  portion  of  the  route  to  the  left  of  0. 

Let  L  be  a  point  such  that  OL  =  X.  If  Aircraft  a  is  between 
0  and  L,  it  is  now  in  conflict  with  Aircraft  b.  However, 
if  it  is  beyond  (to  the  left  of)  L,  it  is  not  now  in  conflict 
with  b  and  could  not  have  been  in  conflict  in  the  past,  as 
can  easily  be  checked  by  moving  each  aircraft  backward.  How¬ 
ever,  as  b  flies  toward  D  and  a  flies  toward  0,  the  distance 
between  them  becomes  shortest  when  the  distance  of  a  from  0 
along  AO  is  equal  to  the  distance  of  b  from  0  along  OD .  (This 
fact  can  be  proved  either  mathematically  or  by  simple  geometric 
trials.)  We  must  make  sure  that  this  shortest  distance  is  ^  X. 
Let  R  be  the  nearest  point  to  0  such  that  Aircraft  a  being  at 
this  point  when  Aircraft  b  is  at  0,  will  come  at  most  no  closer 
than  X  miles  to  b,  as  both  aircraft  continue  their  flight. 

Now,  referring  to  Figure  C-2,  let  a  be  at  M  and  b  at  N ; 
when  the  aircraft  reach  these  points, 

OM  =  ON  (condition  for  shortest  distance) 

But  since  the  aircraft  are  assumed  to  travel  with  the  same 

speed,  the  distance  RM  covered  by  a  should  be  equal  to 

ON — the  distance  covered  by  b  in  the  same  time  period.  Thus, 

RM  =  ON  =  OM 
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FIGURE  C-2  SHORTEST  DISTANCE  BETWEEN  AIRCRAFT  a  AND  b 
By  assumption,  MN  =  X  miles.  Now, 


(OM) 
l/2  (MN) 


O' 


-  sec  — 


2 


OM  =  l/2(MN) 


sec 


a 

2 


o' 

OR  =  MN  s  “  =  X  sec 
2 


Ot 

2 


Therefore,  the  critical  length  to  the  left  of  0  is  equal  to 
X  sec  Qf/2. 

•  Consider  now  the  portion  of  the  route  to  the  right  of 
0.  Let  Q  be  a  point  X  miles  from  0.  If  an  aircraft 
a  is  somewhere  between  0  and  Q,  it  is  now  in  conflict 
with  b.  If  it  is  beyond  Q,  then  it  is  not  now  in  conflict 
with  b  and  cannot  come  into  conflict  in  future.  However, 
considering  the  past  and  applying  a  reasoning  similar  to 
that  used  to  calculate  the  length  OR  on  the  left-hand  side 
for  the  future,  a  length  OP  equal  to  OR  exists  on  the  right 
hand  side  such  that  Aircraft  a  in  the  strip  QP  has  been  in 
conflict  with  Aircraft  b,  as  can  be  checked  by  moving  Air¬ 
craft  a  now  between  Q  and  P,  and  Aircraft  b,  now  at  O,  back 
ward.  Thus  the  critical  length  on  the  right-hand  side  is 
also  X  sec  cy/2.  Denoting  the  length  RP  by  L  ,  the  critical 
length,  we  have: 


L 

c 


2X  sec  “  nmi 
2 
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a .  Expected  Number  of  Potential  Conflicts 

Every  aircraft  of  flow  f  along  AB  will  occupy  the  critical 
Lc  1  ' 

length  L  for  a  period  of  —  hours.  Therefore,  during  a  period  of  1 
c  V 

hour,  the  critical  length  L  will  be  occupied  by  flow  f  for  a  period 
7  c  1 

of  f  —  hours.  The  number  of  aircraft  along  CD  arriving  at  the 

1  V  Lc 

intersection  during  this  period  will  be  f  —  f  .  It  is  not  difficult 
L  1  V  2 

c 

to  see  that  f  f^  is  the  expected  number  of  potential  conflicts  per 

hour.  Denoting  by  E  the  expected  number. of  potential  conflicts  per 

c 

hour,  and  substituting  for  L  ,  we  get 

c 


f 

12  Q* 

E  =  2X  — ■ —  sec  —  conflicts  per  hour, 
c  V  2 


(C-l) 


The  above  expression  indicates  that,  for  fixed  values  of  X,  f  ,  f  ,  and 

x  Z 

V,  the  expected  number  of  conflicts  increases  as  <y  increases. 


b.  Saturation  Flow  Rate 


An  interesting  interpretation  of  the  above-noted  relationship 
is  as  follows: 

Let  f^  ^  f2-  This  entails  no  loss  of  generality.  For  fixed 
values  of  O'  and  X,  the  value  of  f-^  that”  makes  the  expression 
2X  f^/V  sec  Oi/2  equal  to  unity  will  block  the  intersection 
continuously  and  can  be  regarded  as  the  saturation  flow 
rate  at  the  intersection.  With  this  flow  rate  on  one 
route,  every  aircraft  on  the  crossing  route  will  experience 
a  conflict. 

Thus,  the  saturation  flow  rate  at  an  intersection  is: 


f 

s 


V_ 

2X 


Of 


cos  - 


(C-2) 
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c .  Conflict  Intensity  Index 

A  convenient  method  to  characterize  the  conflict  characteristics 

of  an  intersection  would  be  to  normalize  Equation  (C-l)  with  reference 

to  the  product  f  f  .  Thus,  defining  fi/f  f  as  the  conflict  intensity 
12  12 

index,  I,  we  have 

1  =  Ec/flf2  =  f  Se°  f  .  '  (C_3) 

Note  that  I  =  l/f  .  For  any  specified  values  of  X  and  V,  a  normalized 
s 

graph  indicating  the  relationships  between  I  and  O'  and  f  and  o  can 

s 

be  established.  These  curves  can  be  used  as  a  convenient  reference  to 

calculate  the  number  of  predicted  conflicts  when  O',  f  ,  and  f  are 

1  2 

known. 

In  Figure  C-3,  the  relationship  between  f  ~  ot  and  I  ~  a  is 

s 

shown  for  X  =  5  nmi  and  V  =  500  mph. 

2 .  Mathematical  Model  for  Conflict  Duration 

The  potential  conflict  prediction  model  based  on  radar  separation 
rules,  as  described  above,  gives  only  the  number  of  potential  conflicts 
to  be  expected  per  unit  time  but  does  not  give  any  indication  about  the 
expected  duration  of  these  conflicts.  The  purpose  of  the  discussion  be¬ 
low  is  to  develop  an  analytical  expression  giving  the  expected  duration 
of  conflicts  as  a  function  of  the  intersecting  angle  ck  and  the  rates  of 
flow  along  the  intersecting  routes...  It  is  shown  that,  for  the  lower 
values  of  O',  the  number  of  conflicts  is  less  but  the  expected  duration 
of  the  conflicts  is  longer,  compared  with  the  higher  values  of  o  f  for 
which  the  number  of  conflicts  is  greater  but  the  expected  duration  of  the 
conflicts  is  shorter.  Some  relationships  between  conflict  duration  and 
intersection  capacity  are  also  developed. 
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FIGURE  C-3  CONFLICT  INTENSITY  INDEX  AND  SATURATION  FLOW  RATE  AS  A 
FUNCTION  OF  a,  ASSUMING  X  =  5  MILES  AND  V  =  500  mph 


Expected  Time  in  Conflict 


Consider  Routes  AB  and  CD  intersecting  at  an  angle  O',  as 
shown  in  Figure  C-4 .  Suppose  Aircraft  b  is  now  at  O  traveling  toward  D. 


FIGURE  C-4  TWO  INTERSECTING  ROUTES  FOR  EVALUATION  OF  CONFLICT  DURATION 
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Let 


OR  =  OP  =  X  sec  a/2, 

Where  X  =  the  radar  separation  minimum. 

We  know  (from  Section  2)  that  if  Aircraft  a  along  AB  is  anywhere  in  RP , 
then  one  of  the  following  is  true  in  relation  to  Aircraft  b: 

•  It  is  now  in  conflict 

•  It  has  been  in  conflict 

•  It  will  be  in  conflict. 

If  Aircraft  a  is  outside  the  section  RP,  then  no  conflict  between  Air¬ 
craft  a  and  b  is  occurring,  has  occurred,  or  will  occur.  We  concern 
ourselves  with  the  condition  when  b  is  at  0  and  a  is  somewhere  in  RP . 

Our  assumption  is  that  any  point  along  RP  is  an  equally  likely  one  for 
Aircraft  a  to  be.  Thus,  consider  a  point  H  a  distance  of  y  miles  away 
from  0,  and  assume  that  Aircraft  a  is  at  H  when  Aircraft  b  is  at  O.  The 
duration  of  conflict  in  this  situation  will  be  the  period  starting  from 
the  time  Aircraft  a  and  Aircraft  b  first  come  within  X  nmi  of  each  other 

until  the  time  when  the  distance  ceases  to  be  less  than  or  equal  to 

X  nmi  . 

Let  both  Aircraft  b  and  a  be  moved  back  in  time,  so  that  b 

is  at  S  and  a  is  at  T,  such  that  OS  =  TH  =  Z .  The  largest  value  of  Z 

such  that  the  length  TS  is  just  X  will  be  the  situation  when  a  and  b 

first  entered  into  potential  conflict  situation.  From  then  onward  the 

distance  between  a  and  b  continued  to  be  less  than  X  until  a  reached  the 

point  F  and  b  reached  G,  with  OF  =  Z  and  OG -=  Y  +  Z.  Thus  the  duration 

of  conflict  is  the  time  spent  in  traveling  the  distance  Y+Z+Z=2Z+Y 

2Z  +  Y 

i.  e.,  the  conflict  duration  time  is.  hours. 

V 

We  now  develop  an  expression  for  the  length  2Z  +  Y .  Referring 
to  Figure  C-4 ,  let  Z  be  such  that  ST  =  X  miles;  then, 

2  2  2 
(Y  +  Z)  +  (Z)  -  2(Y  +  Z)  Z  cos  O  =  X  ,  (03) 
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or 


or 


or 


(2-2  cos  a)  Z' 


( 2Y  -  2Y  cos  a)  Z 


(Y2  - 


2 

X  )  =  0 


.  (C-4) 


2 

Z  +  YZ 


(Y 


x2) 


2(1  -  cos  O') 


0 


2 

Z  +  YZ 


2  O' 

4  sin  ” 
2 


0 


a 


cosec 


V  2  2  2  a 

X  -Y  cos  - 


(C-5) 


We  choose  the  plus  sign,  since  a  negative  value  of  Z  would  indicate 
going  forward  in  time,  the  aspect  we  have  already  taken  into  account 
by  adding  a  length  Z  to  the  right  of  0. 
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Thus 


Z  = 


O'  4/  2  2  2  O' 

-  Y  +  cosec  “  W  X  -  Y  cos  T 
2  J  2 


o’  4/  2  2  2  O' 

and  2Z  +  Y  =  conflict  duration  distance  =  cosec  “  w  X  -  Y  cos  ~ 


04/2  2  cv  2 

_  cot  “  W  X  sec  —  -  Y 
2  f  2 


(06) 


The  expected  conflict  duration  distance  D  is  given  by 


D  = 


/ 

2  o 


O' 

X  sec  — 
2 


V 


a  \  2  2  o'  2 


cot  —  ^  X  sec  —  -  Y  dy 


(C-7) 


O' 

cot  — 
_ 2_ 

O' 

X  sec  — 
2 


v 


2  2  a 

X  sec  — 
2 


-  Y 


2  2  a  -1  /  Y  \ 

+  X  sec  —  sin  I -  1 

2  V  f  j 


a 

X  sec  - 

2 


a 

cot  — 
2 


2  X  sec  — 
2 


(n  2 .  2  ^  | 

—  X  sec  —  I  , 
2  2  / 


tt  a 

D  -  7  X  cosec  — 
4  2 


(C-8) 


Assuming  an  average  speed  of  V  mph,  the  expected  conflict  duration  time 
per  conflict  T  is 


TT  X  O' 

T  =  ”  ~  cosec  ~  hours. 


(C-9) 
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b .  Conflict  Duration  Time  per  Hour 

Multiplying  the  expected  number  of  conflicts  by  T  gives  the 
total  time  (as  'a  fraction  of  an  hour)  per  hour  during  which  aircraft' 
along  Routes  AB  and  CD  would  remain  in  conflict  if  nothing  were  done 
by  ATC .  Thus, 

2  f  f 

TT  X  O'  1  2  <y 

H  =  conflict  time  duration  per  hour  =  -  —  cosec  “  —  X  sec  — 


=  flf2 

V 


a  ot 

2  sin  t  cos  t 
2  2 


H  =  tt  —  f  f  cosec  a  .  (C-10) 

u  1 


Let  us  define  conflict  duration  index  h  to  be  =  H/f  f  ,  so 

x  £ 


2 

X 

h  =  tt  —  cosec  &  .  (Oil) 

V 


c .  Some  Typical  Numerical  Results 

Let  X  =  5  nmi  and  V  =  500  mph .  With  these  values,  the  values 
of  T  and  h  are  calculated  for  various  values  of  a.  in  Table  C-l. 

For  the  design  of  a  route  network,  it  appears  that  the  number 
of  potential  conflicts  alone  may  not  be  a  sufficiently  complete  attribute 
The  time  duration  of  the  conflicts  should  also  be  considered.  It  is 
evident  that,  even  though  the  number  of  conflicts  increases  with  o', 
the  average  conflict  duration  time  per  hour  decreases  as  O'  increases. 
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Table  C-l 


AVERAGE  CONFLICT  DURATION  AND  CONFLICT  DURATION 
INDEX  AS  A  FUNCTION  OF  O' 


a 

BH 

■i 

2 

X 

h  =  tt  —  cosec  <y 

2 

V 

5 

0.183 

(ID 

.0036 

10 

.09 

(5.4) 

.0018 

30 

.03 

(1.8) 

.00063 

60 

.016 

(0.94) 

.00036 

90 

.011 

(0.67) 

.00031 

Appendix  D 


AIR  ROUTE  AND  INTERSECTION  CAPACITY  MODELS 

1.  Capacity  of  a  Single  Air  Route  when  no  Other  Route 
Intersects  and  no  Passing  is  Permitted 

Suppose  that  the  speeds  of  aircraft  using  a  common  air  route  vary 

between  V  and  V  ,  V  being  >V  .  Let  the  route  length  be  i.  Let  the 
1  2  1  2 

separation  minimum  be  X  miles;  i.e.,  two  aircraft,  one  behind  the  other, 
should  not  come  closer  than  X  miles  at  any  time.  The  worst  condition 
in  terms  of  one  aircraft  following  another  is  when  an  aircraft  flying 
at  V  mph  is  followed  by  an  aircraft  flying  at  V  .  So  that  the  faster 
aircraft  will  not  overtake  the  slower  aircraft  before  the  slower  has 
reached  the  destination  D,  the  faster  aircraft  must  be  released  from 
the  origin  0  at  a  time  T  hours  after  the  slower  aircraft,  where  T  has 
to  have  a  certain  minimum  value.  We  first  calculate  this  value. 

Referring  to  Figure  D-l,  let  Aircraft  a  be  a  distance  m  miles  away 
from  0  after  T  hours,  when  Aircraft  b  is  released  from  0.  Now 

m  =  TV  .  (D-l) 

2 

- e - - - — 

•  I  —  • 

D  a  0 

-m - m - &*>■ 

SA-1096-20 

FIGURE  D-1  A  SINGLE  AIR  ROUTE  FOR  CALCULATING  OVERTAKING  TIME 
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If  Aircraft  a  reaches  D  when  Aircraft  b  is  X  miles  away  from  D,  no  con¬ 
flict  will  ensue,  since  b  will  already  have  reached  its  destination. 
Thus  we  must  have 


(D-2) 


(D-3) 


T 


U  -  X) 
V 

1 


(D-4) 


4V  -  v  U  -  x) 


V1X2 


(D-5) 


If  more  than  two  aircraft  follow  each  other  during  a  period  of,  say, 

1  hour,  the  worst  sequence  in  terms  of  one  aircraft  overtaking  another 
would  be  as  shown  in  Figure  D-2.  We  denote  by  S  the  slowest  and  by  F 


h— T— h-'mH 


SA-1096-21 

FIGURE  D-2  WORST  SEQUENCE  OF  SLOW  AND  FAST  AIRCRAFT 

the  fastest  aircraft.  Theoretically  a  slow  aircraft  following  a  fast 

aircraft  can  be  released  .from  the  origin  immediately  after  the  release 

of  the  other.  However,  according  to  ATC  rules,  there  has  to  be  a  lapse 

of  a  certain  minimum  time  t  before  the  release  of  the  second  aircraft. 

m 
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Thus,  in  Figure  D-2,  going  from  right  to  left,  the  time  separation 

between  S  and  F  is  t  and  that  between  F  and  S  is  T,  It  can  easily  be 

m 

seen  that  there  is  one  aircraft  in  each  of  the  time  spans  (T/2  4-  t  /2) . 

m 

Therefore,  the  average  number  of  aircraft  that  can  be  released  along 
Route  OD  is 


aircraft  per  hour 


2 


[*V 

t  +  - - 

m 


V2  U  ’  X)] 


V  V 
1  2 


aircraft  per  hour 


(D-6) 


Let  us  denote  this  number  of  aircraft  as  A  ,  the  minimum  capacity  of 

min 


the  air  route: 


A 

min 


2 


aircraft  per  hour 


(D-7) 


Here  A  is  the  minimum  theoretical  capacity,  assuming  the  worst  se- 
min 

quence  of  slow  and  fast  aircraft  when  no  passing  is  allowed.  Actually, 
the  worst  sequence  is  very  unlikely  to  occur,  and  passing  of  aircraft 
is  allowed.  A  more  realistic  measure  of  capacity  is  considered  in  the 
concluding  remarks  on  this  model. 


2.  Some  Typical  Numerical  Results 
Let 

t  =  2  minutes  =  1/30  hour 
m 

X  =  20  miles 
=  300  mph. 
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Then, 


2 


30  V  600 

2 


Values  of  A  for  several  values  of  V  are  shown  for  i  =  500  and  1.000 
min  2 

below.  These  values  are  also  plotted  in  Figure  D-3. 

A  ~  V  (V  =  600  mph) 

min  2  1 


A  . 
min 


2 

&  =  500  miles 

i  =  1,000  miles 

600 

30 

30 

550 

13 

9 

500 

8 

5 

400 

5 

2 

300 

3 

1 

3,  Concluding  Remarks  on  the  Capacity  Model  with  No  Passing  Allowed 

A  study  of  Figure  D-3  indicates  that  the  theoretical  minimum  capacity 

of  a  single  route,  when  no  passing  is  allowed,  decreases  very  rapidly 

even  when  the  speed  of  two  aircraft  differs  only  slightly.  Thus,  with 

both  V  and  V  being  600  nmi,  the  capacity  is  30  aircraft  per  hour. 

1  2 

With  V  being  550  nmi  (i.et,  only  8  percent  less  than  the  fast  speed), 

2 

the  capacity  is  reduced  to  only  13,  less  than  half.  These  numbers  also 
indicate  the  amount  of  workload  on  ATC  to  implement  overtaking  of  air¬ 
craft  if  overtaking  is  allowed.  It  is  easily  inferred  that  slight  varia¬ 
tion  in  the  speeds  of  aircraft  will  significantly  increase  the  ATC  work¬ 
load  associated  with  providing  passing  guidance. 
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FIGURE  D-3  VALUES  OF  A  FOR  SEVERAL  VALUES  OF  V9 

A  more  realistic  measure  of  capacity  would  be  to  consider  the  typi¬ 
cal  mix  of  aircraft,  assuming  suitable  distribution,  and  to  establish  a 
capacity  measure,  possibly  on  the  basis  of  allowing  a  certain  number  of 
passing  aircraft  per  unit  length  of  route  per  unit  time,  as  discussed 
below. 


4.  A  Proposed  Definition  of  the  Average  Capacity 
of  a  Single  Route  with  Passing  Allowed 


Let  a  single  route  be  used  by  various  types  of  aircraft  having 

speeds  V  ,  V  ,  V  ,  ...  mph,  V.  being  greater  than  V.  .  Let  f  be  the 
12  3  i  1+ 1  i 

flow  per  hour  of  the  aircraft  having  speed  V. .  We  define  the  average 

capacity  of  the  route  to  mean  the  sets  of  flows  f  ,  f  ,  f  ,  ...  for 

12  3 

which  the  number  of  overtaking  aircraft  per  unit  length  of  route  per 
unit  time  does  not  exceed  a  specified  value. 
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5. 


Derivation  of  the  Route  Capacity  Model  with  Passing  Allowed 


Let  us  first  consider  only  two  types  of  aircraft,  having  speeds 

and  V  ,  with  V  >V  .  Let  the  corresponding  flows  be  f  and  f  per  hour. 

Consider  a  route  length  of  &  miles.  If  we  consider  flow  £  independently, 

the  average  spacing  between  aircraft  with  speed  V  will  be  V  /f  .  The 

2  2  2 

speed  of  the  fast  aircraft  relative  to  the  slow  aircraft  is  (V  -  V  )  mph. 

1  2 

During  a  period  of  1  hour,  the  string  of  fast-moving  aircraft  will  travel 
a  distance  of  (V^  -  V  )  miles  relative  to  the  string  of  slow-moving  air¬ 
craft.  During  this  period,  each  fast-moving  aircraft  will  overtake 

(V  -  V  )  £  /V  slow-moving  aircraft.  In  the  length  H  there  will  be,  on 
12  2  2 

the  average,  i  f^/V^  fast-moving  aircraft.  Thus,  during  a  period  of  1 
hour,  the  total  number  of  overtakings  G  will  be 


G  = 


1  v2(vi  -  y 

V  V 
1  2 


(D-8) 


If  there  were  other  types  of  aircraft  having  speeds  of  V  ,  V  ,  ...  with 

3  4 


V 


>V  >V  ,  the  total  number  of  overtakings  per  hour  would  be: 

1  ^  J 


it. 

G  =  — 1 
T  V 


f(V-V)  f(V-V)  f(V-V) 
212  313  41  4 

- - -  4-  - -  +  -  + 


if 


f  (V  -  V  )  f  (V  -  V  ) 

3  2  3  4  2  4 

-  +  -  + 


(D-9) 
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For  any  specified  value  of  and  i  (e.g.,  two  overtakings  per  hour 
per  100  miles),  the  capacity  of  the  air  route  is  the  set  of  all  the 
combinations  of  the  values  of  f  ,  f  ,  f^,  ...  that  satisfies  Eq.  (D-9) . 


6.  An  Example 

Let  us  consider  the  simple  case  of  two  types  of  aircraft  with 

V  >  V  .  Letting  i  =  500  miles,  we  have 
12  &  ' 


G  = 


50°  ixt2  (Vl  -  V2) 

V  V 
1  2 


(D-10) 


If  G  equals  one  overtaking  allowed  per  hour,  then 


V  V 

f  =  1  2 

fl  2  500  (V 


y 


gives  the  capacity  of  the  air  route  per  500  miles  in  terms  of  the  product 

of  the  two  types  of  aircraft.  Let  V  =  600;  then  the  values  of  f  f 

1  •  ».  12 

for  various  values  of  are  as  shown  in  the  following  table: 


V 


f  f 
1  2 


&  =  500  miles  i  =  250  miles 


550 

500 

450 


13 

6 

3 


26 

12 

7 


Note:  For  =  Mg,  we  must  use  the  capacity  based  on  the  minimum  allow¬ 

able  separation,  e.g.,  if  the  separation  minimum  is  X,  the  capac¬ 
ity  is  V/X  aircraft  per  hour. 
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7. 


A  Special  Case  of  Capacity  Calculations  with  Passing  Allowed 


If  f  .  f  ,  f  .  . . .  are  known  as  fractions  of  the  total  flow  f, 
l'  2’  3’  ’ 

then  Eqt  (D-9)  can  be  rearranged  in  a  more  convenient  form.  Let  f  .  =  k  f 

~  i  i 

then. 


(D-ll) 


f 


G  /i 

T 


0  (k  ,V 

i  i 


) 


aircraft  per  hour 


(D-12) 


where  0(k.,Vi>  is  equal  to  the  expression  in  braces  in  Eq.  (D-ll). 
Equation  (D-12)  gives  the  capacity  of  a  single  air  route  in  terms  of 
(G^/f),  which  is  the  allowable  average  number  of  overtakings  per  hour 
per  unit  length  of  air  route.  The  speeds  of  various  aircraft  and  the 
fraction  of  the  total  flights  for  each  speed  are  assumed  to  be  known, 
(Note:  It  is  assumed  that  V  >  V.^  without  any  loss  of  generality.) 
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8 .  An  Example  for  the  Special  Case 

Let  us  consider  a  mix  of  three  types  of  aircraft : 

Let  V  =  600  V  =  500  V  =  400 

1  2  3 

Let  k  =  0.6  k  =  0.2  k  =  0.2 

12  3 

Let  =  1/250  i  i.e.;  one  overtaking  per  250  mph . 

Substitution  of  the  assumed  values  in  Eq.  (D-12)  would  give 

f  =  5  aircraft  per  hour 

9 .  Concluding  Remarks  on  Capacity  Model  with  Passings  Allowed 

A  definition  of  average  air  route  capacity  based  on  a  certain  allow¬ 
able  number  of  overtakings  per  hour  per  unit  length  of  air  route  has 
been  presented  above.  If  the  mix  of  aircraft  of  various  speeds  is  not 
known,  the  capacity  is,  a  set  of  various  values  of  flows  f  ,  f  ,  f  ,  ... 
fulfilling  a  certain  constraint  (Eq.  D-9) .  If  the  mix  or  fraction  of 
aircraft  of  various  speeds  is  known,  capacity  can  be  calculated  as  a 
single  number  (Eq.  D-12).  Two  examples  have  been  presented  to  illustrate 
the  concepts.  These  concepts  will  be  pursued  further  and  any  new  ap¬ 
proaches  reported  in  due  course. 

10 .  Capacity  at  an  Intersection 

Referring  to  Figure  D-4,  a  possible  definition  of  an  aircraft  capac¬ 
ity  at  the  intersection  of  two  routes  would  be  the  maximum  allowable 

value  of  the  product  f  f  (where  f  and  f  are  the  flows  of  aircraft 

12  1  2 

per  hour  on  the  two  routes)  under  one  of  the  following  restrictions: 
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(1)  The  expected  number  of  potential  conflicts  per  hour 
does  not  exceed  a  certain  specified  number. 

(2)  The  expected  time  the  aircraft  along  the  two  routes 
spend  in  violation  of  the  separation  minimum  does  - 
not  exceed  a  specified  fraction  of  an  hour. 


FIGURE  D-4  TWO  INTERSECTING  ROUTES  FOR  CAPACITY  EVALUATION 


11 .  Intersection  Capacity  Model  under  Restriction  One 

From  Appendix  C,  we  know  that 

E  =  expected  number  of  potential  conflicts  per  hour 
c 

2f  f 

12  Qf 

=  -  X  sec  -  .  (D-13) 

V  2 


Rearranging,  we  get 


f  f 
1  2 


E  V 
c 


2  X  sec  — 
£ 


(D-14) 


Letting  E 

c 


1,  i.e.,  one  conflict  per  hour  allowed,  we  have 


^  ^  -V  O' 

f  f  =  —  cos  — 

1  2  2X  2 


(D-15) 
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Equation  (D-15)  gives  the  capacity  at  an  intersection  as  a  function  of 
average  aircraft  speed  V,  separation  minimum  X,  and  intersection  angle, 
O',  under  the  restriction  that  one  conflict  per  hour  is  acceptable  on  the 
average.  The  capacity  is  directly  proportional  to  the  allowable  number 
of  conflicts  per  hour. 

Assuming  typical  values  for  V  and  X  to  be  600  and  5,  respectively, 
Table  D-l  gives  values  for  capacity  as  a  function  of  intersection  angle 
a.  These  values  are  plotted  in  Figure  D-5  for  convenience. 

Table  D-l 

INTERSECTION  CAPACITY  AS  A  FUNCTION  OF  INTERSECTION  ANGLE 
(One  Conflict  per  Hour  Acceptable) 


Intersection  Angle  o' 
(degrees) 

Capacity  f-^fg 
(unit  conflict  per  hour) 

10 

59 

20 

59 

30 

58 

45 

55 

60 

51 

90 

42 

120 

30 

12,  Intersection  Capacity  Model  under  Restriction  Two 


From  Appendix  C  we  have 


ttX 

H  =  conflict  duration  time  per  hour  =  — -  f  f  cosec  rv  (D-16) 

2  12 


^  ^  HV 

f  f  =  -  sin  ca 

12  2 
ttX 
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FIGURE  D-5  IfSITERSECTIOM  CAPACITY  AS  A  FUNCTION  OF  a 


Letting  H  =  0.01,  i.e.,  the  conflict  duration  time  per  hour  is  1  percent 
per  unit  time,  we  get 

v2 

f  f  =  0.01  —  sin  v  .  (D-17) 

i  2.  £ 

ttX 

In  Table  D-2,  values  of  f  f  as  a  function  of  <y  are  tabulated  for 

1  2 

V  =600  and  X  =  5.  The  results  are  also  plotted  in  Figure  D-5. 

13 .  Concluding  Remarks  on  Intersection  Capacity 

A  study  of  Figure  D-5  indicates  the  interesting  fact  that,  for 
a  <  90°,  the  intersection  capacity  based  on  conflict  number  decreases 
slowly  as  a  increases,  whereas  the  capacity  based  on  conflict  duration 
time  increases  as  <y  increases.  For  c y  >  90°,  both  capacities  decrease 
as  o'  increases.  Further  study  appears  to  be  necessary  to  establish 
which  capacity  definition  corresponds  best  to. real  world  situations. 
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Table  D-2 


INTERSECTION  CAPACITY  AS  A  FUNCTION  OF  INTERSECTION  ANGLE 
(1-percent  Time  Allowed  in  Conflict) 


Intersection  Angle  a 
( degrees) 

Capacity  1^2 

(for  1-percent  conflict  duration) 

10 

8 

30 

23 

45 

32 

60 

39 

90 

45 

120 

39 

150 

23 
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Appendix  E 


ASSIGNING  AN  ALTITUDE  OTHER  THAN  THE  DESIRED  ONE 

The  purpose  of  this  appendix  is  to  show  that,  within  the  framework 
of  the  hemispherical  rule,  altitude  assignment  is  not  a  practicable 
design  variable. 

1 .  Consequence  of  Assigning  Altitudes 
Other  than  that  Desired  by  the  User 

Consider  a  Route  AB  such  that  it  intersects  Routes  A^B  ,  A^B^.,  •••> 

A  B  at  the  common  altitude  level  (e.g.,  31,000  ft),  as  shown  in  Figure 
n  n 

E-l. 


FIGURE  E-l  SEVERAL  INTERSECTING  ROUTES  FOR  ALTITUDE  ASSIGNMENT 


A  possible  method  of  removing  the  intersection  is  to  assign  to 

Route  AB  an  altitude  other  than  that  used  by  Routes  A  B  ,  A  B  ,  .  .  .  , 

1  1 9  2  2’  ’ 

and  so  on.  Frequently,  however,  it  is  desirable  for  a  route  such  as 
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>SJV- 

AB  to  have  the  same  altitude  as  that  of  A^B^,  A^B^,  ’  '  '  ’  on  bas^s 
of  optimum  fuel  consumption  rate  and  so  forth.  Thus,  if  flights  along 
AB  were  assigned  an  altitude  one  level  lower  or  higher,  the  fuel  costs 
would  be  increased. 


a. 


Fuel  Costs 


Let  the  increase  in  fuel  cost  be  $K^  per  100  miles  if  an  alti 
tude  other  than  that  desired  is  used.  If  the  length  of  Route  AB  is  Hy 
and  the  flow  is  Y  aircraft  per  hour,  the  extra  fuel  cost  for  using 
a  nonoptimum,  flight  level  will  be: 


C  =  K  Y  - — -  dollars  per  hour. 
1  1  100 


(E-l) 


b. 


Conflict  Costs 


If  the  desired  altitude  is  used,  this  being  the  same  as  that 

used  by  A  B  ,  A  B  ,  ....  then  potential  conflicts  will  result.  If 

11  2  2  7 

x  ,  x  ,  x  ,  ...  is  the  flow  of  aircraft  per  hour  along  A  B  ,  A  B  ,  ... 
12  3  112  2 

the  number  of  expected  potential  conflicts  per  hour  E  is 

c 


-  (x  -f-  x  +  x  + 
V  1  2  3 


O' 


)  X  sec  — 
2 


(E-2) 


where 

V  =  average  aircraft  speed  (miles  per  hour) 

X  =  separation  minimum 

O'  =  average  angle  of  intersection. 

Assuming  V  =  500,  X  =  5,  and  s  01 /2  ~  1,  we  get 


Y(X1  +  X2  +  V 

50 


(E-3) 


100 


the  cost 


Let  K  be  the  cost  per  conflict  in  dollars,  then  C  , 
2  2 

of  conflicts  per  hour,  is 


V 

50 


(X1  +  X2  +  X3  + 


(E-4) 


Considering 
advantageous 
it  is  not--i 


costs  C  and  C  ,  it  is  clear  that,  if  C  <  C  ,  then  it  is 
12'’  ’  1  2’ 

to  assign  an  altitude  other  than  that  desired;  otherwise 
.  e.,  for  using  undesired  altitudes  the  condition  is 


l  K2  ■ 

V  155  <  To  Y  (X1  +  X2  +  V 


(E-5) 


Cancelling  Y  on  both  sides  and  rearranging, 


K1  2(X!+X2+ 


(E-6) 


or 


^  twice  the  sum  of  flights  per 

1  hour  on  intersecting  routes 

—  <  - 2 -  .  (E-7) 

K  length  of  candidate  route 

for  altitude  modification 


c .  An  Example 

We  now  consider  a  typical  case.  The  direct  route  from  New 
York  to  Los  Angeles  intersects  13  routes  in  the  first  1,000  miles  of 
its  length  at  the  altitude  of  35,000  ft.  The  sum  of  the  flights  on  the 
13  routes  at  35,000  ft  is  typically  50  flights  during  5  hours  of  peak 
traffic,  assuming  that  the  peak  traffic  occurs  in  the  same  hours  (actu¬ 
ally,  the  sum  will  be  less  than  50  because  of  diversity  in  peaking). 
This  gives  about  ten  flights  per  hour  on  the  intersecting  routes.  If, 
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on  the  basis  of  these  figures,  the  route  from  New  York  to  Los  Angeles 
is  to  be  assigned  an  altitude  of,  say,  39,000  ft  for  the  first  1,000 
miles,  thereby  removing  all  intersections  (there  is  no  cross  traffic 
at  this  level),  we  must  have 

K 

_1  2  X  10 

K  1000 

2 

1 

<  50 

i .  e, , 

Cost  per  100  nmi  of  the  altitude  change  1 
Cost  per  conflict  50 


It  has  been  estimated  that  for  a  jet  aircraft,  the  cost  per  100  nmi  of 
an  altitude  change  is  about  $10.  Thus,  under  present  traffic  conditions, 
unless  the  cost  of  a  conflict  were  $500  or  more,  it  would  not  be  profit¬ 
able  to  change  the  altitude  of  the  New  York  to  Los  Angeles  route  from 
35,000  to  39,000  ft  (for  the  first  1,000  nmi).  Similar  remarks  are 
obviously  applicable  to  other  routes. 


2.  Consideration  of  Future  Traffic 


If  the  traffic  increases  by  a  factor  of  n. 


the  number  of  conflicts 


increases  by 
traffic,  the 


a  factor  of  n 
relationships 


K  nY 
1 


t 

100 


< 


.  Therefore,  assuming  an  n-fold  increase  in 
(E-5)  and  (E-6)  will  be  modified  to 


2 

K  n  Y 
2 

-  (x  +  x  +  x  +  ...)  (E-8) 

50  1  2  2 


or 


K1  2n(xi  +  x2  +  X3 

- * - 


(E-9) 
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Considering  the  route  from  New  York  to  Los  Angeles  again,  we  have 


K  <  50 
2 

Letting  K  =  $10  and  K  =  $50  per  conflict,  we  get 
1  2 

10  n_ 

50  50  ’ 

or  n  would  have  to  be  10  or  greater  before  it  would  be  profitable  to 
change  the  altitude  of  the  New  York  to  Los  Angeles  route  from  35,000 
to  39,000  ft.  Similar  remarks  are  applicable  to  other  routes. 

The  justification  of  assigning  an  altitude  other  than  that  desired 
is  a  saving  in  costs.  The  calculations  for  a  typical  route  indicate 
that,  with  assumed  costs  of 

$10  per  100  miles  of  flight  at  an  undesired  altitude,  and 
$50  per  conflict, 

the  present  traffic  intensitiy  would  have  to  increase  at  least  tenfold 
for  altitude  change  to  be  profitable.  Thus,  the  principle  may  be  appli¬ 
cable  to  traffic  conditions  some  time  in  the  far  future,  but  it  does 
not  appear  to  be  of  much  advantage  now,  unless  the  conflict  cost  increases 
significantly. 
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Appendix  F 


THE  COORDINATES  AND  INTERSECTION  ANGLE 
OF  ROUTE  INTERSECTING  POINTS 


I.  Necessary  and  Sufficient  Conditions 

for  Intersection  of  Two  Route  Segments 

a.  Conditions  in  a  Plane 


Let  two  route  segments,  i  and  j,  in  a  plane  be  defined  by  their 

end  points  (x  ,  y.-),  (x  ,  y.0)  and  (x  ,  y_),  (x.  .  y .  n)  ,  as  shown 
ll  ll  7  i2'  i2  jl  jl  j 2  j 2 

in  Figure  F-l. 


FIGURE  F-1  TWO  ROUTE  SEGMENTS  IN  A  PLANE 


Let  a  general  point  on  these  segments  be  expressed  as 


x  =  X  x  +  (1  -  \  ) x 

n  ij  il  ij  12 


y  =  .y. 1  +  (l  -  .)y.0 

n  ij  il  ij  i2 


X  =  U  X-  +  ( I  -  U  )  X 

m  ij  jl  ij  j 2 


y  =  n.  .y..  +  (i  -  m* .  ,)y  0 

.  m  iJ  Jl  ij  j2 


(F-l) 


(F-2) 
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where  \  and  u.  are  two  independent  constants.  If,  for  some  values 
ij  ij 

of  \  and  \i  such  that 

ij  ij  .  •  • 


0  <  X  <1  and  0  <  u.  .  <  1 
ij  U 


we  can  obtain: 


x  =  x  and  y  -  y 
n  m  n  m 


then  this  common  point  will  be 
segments,  If  no  such  values  of 
not  intersect. 


the 

ij 


intersecting  point  of  the  two  route 

and  \i  exist,  then  the  routes  do 
ij 


Now,  solving  for  \  and  \i.  by  using  Cramer’s  Rule,  we  have  : 

ij  ij 


'ij 


(XJ2- 

V 

(XJ2- 

V 

• 

(xu  - 

Xi2) 

<Xj  2 

Xi2> 

(yj2 

V 

(yj2- 

»  ^  ■ 

IJ 

- 

V 

<yj2  ' 

yi2) 

(Xil  - 

(Xj2 

V 

(Xil- 

Xi2> 

(v- 

V 

(*il  - 

(yj  2 

V 

(yil  - 

yi2> 

(yj2  - 

V 

(F-3) 


If  the  solution  of  (F-3)  results  in: 


0  <  \  <  1 
ij 


0  <  ia  <  1 

!j 


then  the  route  segments  intersect;  otherwise,  they  do  not.  This  provides 
us  with  a  necessary  and  sufficient  condition  for  route  intersections  in 
a  plane.  A  corresponding  relationship  using  spherical  coordinates  is 
established  below. 


108 


b.  Conversion  to  Sperical  Coordinates 

Let  the  end  points  of  the  two  route  sections,  i  and  j,  be  given 
on  the  earth’s  surface  in  terms  of  latitude  and  longitude  as  follows: 


Route  i:  0. 


il 


latitude  at  one  end  of  Route  i 


Route  j 


0  =  longitude  at  one  end  of  Route  i 

11 

0  =  latitude  at  the  other  end  of  Route  i 

1  ^ 

0  =  longitude  at  the  other  end  of  Route  i 

12 

0  ,  0.  ,  0  ,  0  correspondingly, 

j  1  jl  j2  j  2 


These  coordinates  can  be  transformed  to  plane  coordinates  by  using 
gnomonic  projection,  which  transforms  great-circle  routes  to  straight 
lines.  Also,  since  the  condition  of  intersection  is  independent  of  the 
coordinate  system,  the  Equation  (F3)  can  be  used  to  establish  whether 
the  two  routes  intersect.  The  gnomonic  transformation  gives  the  fol¬ 
lowing  relationships: 


X  „ 

=  cot 

0., 

cos  9  n 

il 

il 

il 

y 

=  cot 

0 

sin  -  9 

ll 

il 

il 

X 

=  cot 

0  r. 

cos  9  _ 

i2 

i2 

i2 

y 

=  cot 

0. , 

sin  0 

i2 

i2 

i2 

for  Route  i  (F-4) 


J 


and  x.^,  y x  2*  ^  2  corresP°nc^n£^y •  We  can  use  these  values  in 

Equation  (F3)  to  calculate  X  and  |i  and  thus  establish  whether  Routes 

ij  ij 

i  and  j  intersect . 


2. 


Coordinates  of  the  Intersection  Point 


If  the  two  routes  intersect,  the  latitude  0  and  longitude  0  of 

ij  i  j  • 

the  intersection  point  can  be  calculated  as  follows: 
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Let 


x  =  cot  0  .  cos  9  .  . 
ij  iJ  ^ 


v  =  cot  0  sin  9  .  . 
l  j  ij 


be  the  gnomonic  projections  of  the  intersection  point.  From  Equation 
(F“l)  we  have 


=  x  =  x  =  cot  0  cos  9  =  A-  x_+(l-A..)x 

ij  n  m  ij  ij  ij  il  ij 


v  =  v  =  V  =  cot  0  sin  9  =  A  y  +  ( 1 -A  )  y 

ij  n  m  ij  ij  ij  il  ij  i- 


,  (F-5) 


which  gives 

tan  6 

ij 


0  =  tan 

ij 


0  =  cot 

ij 


A 

ij 

yu 

+ 

(1-A.  .) 
ij 

yi2 

A 

x  n 

+ 

(1-A  .  .) 

X 

ij 

il 

ij 

i2 

A  * 

y 

+ 

(1-A  .) 

y.  0 

ij 

il 

IJ 

i2 

A 

X . 

+ 

(1-A.  ) 

x  . 

ij 

ll 

ij 

i2 

A 

y . 

+ 

(1-A  ) 

y.n 

ij 

il 

1  J 

i2 

sin  9 


ij 


(F-6) 


(F-7) 


3 .  Angle  of  Intersection 

Referring  to  Figure  F-2,  let  C  be  the  center  of  a  sphere  representing 

the  earth.  Consider  two  intersecting  routes'  i  and  j,  with  0,0  as 

ij  ij 

the  latitude  and  longitude  of  the  intersection  point,  as  calculated  above. 
The  angle  O'  between  the  two  intersecting  routes  is  the  same  as  the 

1  j 

angle  between  the  normals  to  the  two  planes,  one  containing  the  points 

(0^,  9^),  (0.j,  9_),  and  C  ancJ  other  containing  the  points  .(0.^, 

9..),  (0.  .),  and  C. 

Jl  iJ  iJ 
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FIGURE  F-2  INTERSECTION  ANGLE  BETWEEN  TWO  GREAT-CIRCLE  ROUTES 


Since  the  angle  between  the  planes  under  consideration  is  not  dependent 
upon  the  magnitude  of  the  radius  of  the  sphere,  let.  the  radius  be  as¬ 
sumed  to  be  unity  for  convenience.  Define  a  tangent  plane  AOB  as  the 
Z  =  0  plane  on  which  gnomonic  projections  will  be  considered.  Let  the 
tangent  point  0  be  the  origin,  so  that  the  coordinates  of  C  are  (0,0,1). 

For  the  gnomonic  projections,  the  two  planes  mentioned  above  can 
be  expressed  in  determinant  form  as  follows: 


Plane  1  containing  (0  ,0  ),  (0  ,9  ),  and  C 

il  il  ij  ij 


x  x  x  .  0 

il  ij 


y  y,,  y.  .  o 

ii  ij 

z  o  o  o 


l 


=  o, 


(F-8) 


where  x  and  y  are  given  by  (F-4),  and  x  and  y  are  given 
il  il  ij  ij 

by  (F-5). 


Ill 


Plane  2  containing  (0.,,  9  ),  -(0. 

Jl  Jl  iJ 


9  ) ,  and  C  : 

ij 


X 

xjl 

X 

ij 

0 

y 

yjl 

y.  . 
ij- 

0 

z 

0 

'  0 

0 

1 

1 

i 

1 

The  normals  to  the  above  planes  have  components  proportional  to 
the  cofactors  of  x,  y,  and  z  in  the  above  determinants.  Thus,  the 
normal  vectors  are  given  by: 


(F-10) 


Now,  by  using  standard  relationships  from  vector  analysis,  the  angle 
o'  between  these  two  vectors  X  and  X  is  given  by : 

ij  i  j 
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a 

ij 


which  expands  to : 


-1 


O'  .  .  =  cos 
IJ 


V 


(y .  .  -  y . , ) (y .  .  -  y )  +  (x .  .  -  x  ) (x .  .  -  x  ) 
ij  ii  ij  ji _ ij  ii  ij  ji 

2  2  2 
{(y.  .  -  y.,)  +  (x  -  x  )  +  (x  y..  -  y  x  )  } 

ij  ii  ij  ii  ii  ij  ii  ij 


+  (x  y.  .  -  y  ,x-  )(x  y  -  y.  x.  .) 

ii  ij  il  ij  jl  ij  jl  ij _ 

"  2  •  2  2 
l(y, ,  -  y„)  +  (X^  J  -  x_)  +  (x^,  y,  j  -  y.,x;  J) 


ij  jl 


ij  jl 


jl  ij  Jl  ij 


Since  x^ ,  y^  ...  have  already  been  computed,  the  above  expression 
easy  to  calculate. 


(F-12) 


(F-13) 


} 


is 


113 


Appendix  G 


ESTABLISHING  THE  TRAFFIC  DEMAND  MATRIX  A 


The  purpose  of  this  appendix  is  to  describe  how  the  traffic  demand 
matrix  A,  needed  as  an  input  to  the  computer  program,  can  be  established 
from  the  existing  Peak  IFR  data  tapes. 

The  Peak  IFR  demand  data  tapes  have  been  arranged  in  various 
formats  by  NAFEC .  Two  of  the  tapes,  suitable  for  establishing  the 
demand  matrix  A,  are  the  following: 

(1)  The  tape  containing  a  listing  of  the  number  of  aircraft 
per  day  requesting  various  altitudes  between  various 
communities . 

(2)  The  tape  containing  the  count  of  aircraft  between  various 
communities  for  each  hour  during  24  hours  (without  any 
reference  to  the  altitude). 

From  the  first  tape,  the  number  of  aircraft  desiring  various  altitudes 
can  be  computed  as  fractions  of  total  traffic  flow  between  the  respective 
communities.  For  example,  let  the  number  of  aircraft  desiring  altitudes 
of  28,000,  31,000,  and  35,000  ft  between  communities  i  and  j  be  2,  4, 
and  4  respectively.  If  these  three  are  the  only  altitudes  under 
consideration,  then: 

A  fraction  '  =  0.2  of  the  total  flow  of  10  a/c 

2+4+4  1 

desires  the  altitude  28,000  ft, 

0.4  desires  the  altitude  31,000  ft, 
and  0.4  desires  the  altitude  35,000  ft. 

From  the  second  tape,  the  highest  hourly  count  of  traffic  between 
i  and  j  can  be  easily  found.  If  this  highest  count  is  multiplied  by 
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the  various  fractions  calculated  above,  we  will  obtain  the  number  of 
aircraft  desiring  various  altitudes  during  the  peak  hour  traffic  activity. 
Thus,  if  the  highest  hourly  count  between  i  and  j  is  5,  then  it  can  be 
inferred  that : 

1  aircraft  desires  the  altitude  28,000  ft, 

2  aircraft  desire  the  altitude  31,000  ft, 

2  aircraft  desire  the  altitude  35,000  ft. 

Similar  calculations  can  be  performed  for  other  community  pairs  and  the 
demand  matrix  A  established  for  the  desired  routes. 

Study  of  the  actual  data  indicated  that  the  hourly  counts  are 
typically  0,  1,  or  2  aircraft  per  hour.  Multiplying  the  highest  count 
of  2  with  various  fractions  would  frequently  give  the  aircraft  flow  per 
hour  at  various  altitudes  in  fractions,  i.e.,  the  demand  matrix  A  will 
contain  fractions.  There  is  nothing  wrong  with  fractional  aircraft 
flows  as  far  as  the  computation  of  network  attributes  is  concerned. 
However,  to  make  the  demand  matrix  contain  only  integer  numbers,  a 
convenient  method  is  to  multiply  the  fractions  with  a  suitable  factor 
h,  e.g.,  5  and  round  off  the  numbers.  The  numbers  thus  resulting  can 
be  interpreted  as  the  demand  during  h  hours.  In  the  preliminary  test 
runs,  a  value  of  5  was  used  for  h. 
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Appendix  H 


A  PRELIMINARY  LIST  OF  COMMUNITY  PAIRS  EXCHANGING  10  OR  MORE 
IFR  FLIGHTS  PER  DAY,  USING  ALTITUDES  OF  18,000  FEET  OR  ABOVE 

The  list  of  community  pairs  that  starts  on  the  next  page  was  prepared 
on  the  basis  of  the  1969  Peak  IFR  data  supplied  to  SRI  by  NAFEC .  The 
"Enroute  IFR  Air  Traffic  Survey"  peak  day  fiscal  year  1968,  published  by 
the  FAA's  Office  of  Management  Services  was  also  used  as  a  guide. 

The  list  has  been  arranged  so  that  the  reference  communities  are 
always  to  the  east  of  the  paired  communities.  Only  the  pairs  within 
US48  are  considered. 

The  paired  communities  are  listed  in  two  parts.  In  Part  1  are 
those  communities  that  are  more  than  400  miles  apart  and  whose  traffic 
is  expected  to  employ  basically  altitudes  of  28,000  ft  or  above.  In 
Part  2  are  listed  those  communities  that  are  less  than  400  miles  apart 
and  are  expected  to  employ  essentially  altitudes  between  18,000  ft  and 
28,000  ft.  Both  parts  thus  apply  only  to  flights  at  or  above  18,000  ft. 
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Part  1.  Pairs  employing  essentially  altitudes  of  28,OOOf  or  above 
(Exchange  of  10  or  more  flights  per  day) 


Serial 

Number 

Reference 

Community 

Designation 

Paired  Community 

Distance 

Flight 

Designation 

(App.  Statute  Miles) 

Exchange/Day 

1' 

JFK 

MDW 

713 

194 

2 

MIA 

1092 

94 

3 

DET 

482 

66 

4 

LAX 

2451 

55 

5 

SFO 

2571 

47 

6 

CLE 

403 

38 

7 

ATL 

748 

35 

8 

STL 

8  75 

34 

9 

DAL 

13  74 

23 

10 

MSP 

1018 

19 

11 

HOU 

1420 

17 

12 

CMH 

570 

16 

13 

NEW 

1171 

16 

14 

PBI 

1030 

13 

15 

TPA 

1005 

13 

16 

DEN 

1631 

12 

17 

MEM 

962 

12 

18 

SDF 

660 

12 

19 

CLT 

533 

11 

20 

IND 

646 

11 

21 

SEA 

2400 

11 

22 

MDW 

LAX 

1745 

80 

23 

SFO 

1858 

69 

24 

DEN 

920 

60 
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Part  1  -  (Continued) 


Serial 

Number 

Reference 

Comm uni ty 
Designation 

Designation 

Paired  Community 

Distance 

(App«  Statute  Miles) 

Flight 

Exchange /Day 

25 

MDW 

MKC 

414 

60 

26 

DAL 

803 

38 

27 

SEA 

1737 

28 

28 

LAS 

1524 

20 

29 

OMA 

432 

20 

30 

PHX 

1453 

17 

31 

PDX 

1758 

14 

32 

HOU 

937 

12 

33 

TUL 

598 

12 

34 

MEM 

482 

10 

35 

DET 

LAX 

1983 

17 

36 

MSP 

543 

13 

37 

SFO 

2090 

13 

38 

ATL 

596 

10 

39 

HOU 

LAX 

1370 

16 

40 

MIA 

MDW 

1188 

37 

41 

ATL 

604 

36 

42 

NEW 

669 

17 

43 

CLE 

1087 

13 

44 

DET 

1152 

13 

45 

PIT 

1010 

13 

46 

STL 

1070 

10 

Part  1  -  (Continued) 


Serial 

Number 

Reference 

Community 
Des igna tion 

Designation 

Paired  Community 

Distance 

(App.  Statute  Miles) 

Flight 

Exchange/Day 

4  7 

CLE 

LAX 

2049 

14 

4  3 

STL 

492 

13 

49 

ATL 

554 

10 

50 

BOS 

MDW 

851 

45 

51  ' 

MIA 

1255 

22 

52 

LAX 

2596 

18 

53 

DET 

613 

16 

54 

CLE 

551 

15 

55 

PIT 

494 

14 

56 

SFO 

2699 

10 

57 

SEA 

SFO 

678 

44 

58 

DCA 

MDW 

597 

43 

59 

ATL 

543 

20 

60 

MIA 

923 

18 

61 

LAX 

2300 

17 

62 

STL 

712 

17 

63 

DAL 

1185 

13 

64 

BNA 

450 

12 

65 

MSP  * 

930 

10 

66 

TYS 

430 

10 
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Part  1  - 

(Continued) 

Serial 

Number 

Reference 

Paired  Community 

Flight 

Exchange /Day 

Community 

Designation 

Designation 

Distance 

(App.  Statute  Miles) 

67 

DEN 

LAX 

831 

38 

68 

SFO 

949 

24 

69 

SEA 

1020 

16 

70 

PHX 

586 

16 

71 

LAS 

606 

12 

72 

PDX 

985 

12 

73 

DAL 

LAX 

1240 

38 

74 

SFO 

1483 

24 

75 

ELP 

572 

18 

76 

DEN 

663 

12 

78 

PHX 

LAS 

418 

30 

79 

SFO 

653 

20 

80 

PIT 

MDW 

410 

38 

81 

ATL 

521 

17 

82 

PHL 

MDW 

666 

40 

83 

DET 

443 

23 

84 

MIA 

1019 

23 

85 

ATL 

666 

22 

86 

LAX 

2394 

18 

87 

STL 

LAX 

1589 

24 

88 

DEN 

796 

14 

125 


Part  1  -  (Continued) 


Serial 

Number 

Reference 

Community 

Designation 

Designation 

Paired  Community 

Distance 

(App,  Statute  Miles) 

Flight 

Exchange /Day 

89 

STL 

DAL 

547 

14 

90 

SEA 

1705 

10 

91 

LAX 

SEA 

959 

26 

92 

ATL 

MDW 

587 

40 

93 

NEW 

424 

29 

94 

STL 

467 

22 

95 

DAL 

721 

22 

96 

HOU 

701 

17 

97 

LAS 

SF0 

417 

32 

98 

TPA 

ATL 

416 

30 

99 

MDW 

1003 

17 

100 

BAL 

MDW 

606 

23 

101 

ATL 

577 

14 

102 

MIA 

955 

12 

103 

LAX 

2320 

14 

104 

MKC 

DEN 

558 

25 

105 

LAX 

1356 

16 

106 

DAL 

451 

12 

107 

BUF 

MDW 

454 

21 
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Part  1  -  (Continued) 


Serial 

Number 

Reference 

Community 

Designation 

Designation 

Paired  Community 

Distance 

(App.  Statute  Miles) 

Flight 

Exchange /Day 

100 

ABQ 

LAS 

484 

10 

109 

LAX 

664 

14 

110 

SLC 

LAX 

590 

16 

111 

LAS 

508 

16 

112 

SFO 

600 

16 

113 

MEM 

DAL 

420 

22 

114 

HOU 

478 

12 

115 

NEW 

DAL 

443 

23 

116 

TUS 

LAX 

450 

14 

117 

CLT 

MDW 

584 

10 

118 

MSP 

DEN 

700 

14 

119 

LAX 

1530 

10 

120 

MKC 

413 

12 

121 

OMA 

DEN 

488 

15 

122 

OGD 

SAC 

535 

12 

123 

HAR 

SAN 

910 

10 
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Part  1  -  (Concluded) 


Serial 

Number 

Reference 

Community 

Designation 

Designation 

Paired  Community 

Distance 

(App.  Statute  Miles) 

Flight 

Exchange/Day 

124 

ICT 

DEN 

437 

13 

125 

ROC 

MDW 

520 

14 

Total  exchange  of  flights  per  day 
between  community  pairs  of  Part  1  =  2,906 
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Part  2.  Pairs  employing  essentially  altitudes 
between  18,000  ft.  and  28,000  ft. 

(Exchange  of  10  or  more  flights/day) 


Serial 

Number 

Reference 

Community 

Designation 

Designation 

Paired  Community 
Distance 

(App.  Statute  Miles) 

Flight 

Exchange/Day 

1 

JFK 

PIT 

3l7 

53 

2 

DCA 

205 

46 

3 

BUF 

292 

31 

4 

ROC 

250 

27 

5 

SYR 

194 

11 

6 

MDW 

MSP 

355 

84 

7 

STL 

262 

63 

8 

DSM 

309 

19 

9 

CID 

199 

11 

10 

DET 

MDW 

238 

80 

11 

MKE 

252 

17 

12 

HOU 

DAL 

225 

53 

13 

NEW 

316 

44 

14 

SAT 

189 

14 

15 

MIA 

TPA 

205 

54 

16 

CLE 

MDW 

308 

50 

17 

BOS 

PHL 

271 

31 

18 

JFK 

188 

19 
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Part  2  -  (Continued) 


Serial 

Number 

Reference 

Community 

Designation 

Designation 

Paired  Community 

Distance 

(App.  Statute  Miles) 

19 

BOS 

BAL 

360 

20 

SYR 

264 

21 

DCA 

398 

22 

DCA 

CLE 

306 

23 

DET 

396 

24 

PIT 

192 

25 

DAY 

380 

26 

CMH 

329 

27 

DEN 

SLC 

371 

28 

ABQ 

334 

29 

DAL 

SAT 

252 

30  ■ 

OKC 

190 

31 

LBB 

292 

32 

PHX 

LAX 

357 

33 

SAN 

299 

34 

ONT 

323 

35 

PIT 

DAY 

227 

36 

DET 

205 

37 

CVG 

257 

Flight 

Exchange/Day 

17 

10 

56 

21 

20 

20 

16 

10 

23 
16 

24 
12 
23 

39 

28 

20 

13 

13 

12 
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Part  2  -  (Continued) 


Seri al 

Number 

Reference 

Community 
Des ignation 

Designation 

Paired  Community 

Distance 

(App .  Statute  Miles) 

Flight 

Exchange/Day 

38 

PHL 

PIT 

259 

32 

39 

CLE 

360 

21 

40 

STL 

MKC 

238 

35 

41 

TUL 

361 

14 

42 

LAX 

SAC 

361 

34 

43 

PDX 

357 

18 

44 

> 

FAT 

204 

17 

45 

SJC 

307 

70 

46 

SFO 

337 

237 

47 

ATL 

SDF 

319 

21 

48 

MEM 

337 

18 

49 

BNA 

214 

17 

50 

CVG 

MDW 

252 

29 

51 

STL 

309 

12 

52 

EKE 

MSP 

299 

28 

53 

DAY 

MDW 

232 

27 

54 

IND 

STL 

262 

25 

55 

GEG 

SEA 

229 

24 

56 

PDX 

310 

14 
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Part  2 


(Continued) 


Serial 

Number 

Reference 

Community 
Des ignation 

Designation 

Paired  Community 

Distance 

(App*  Statute  Miles) 

Flight 

Exchange /Day 

57 

MKC 

ICT 

177 

13 

58 

BUF 

DET 

216 

17 

59 

PIT 

178 

11 

60 

ABQ 

PHX 

330 

20 

61 

SLC 

BO  I 

296 

19 

62 

RNO 

SFO 

185 

20 

63 

CMH 

MDW 

252 

19 

64 

MEM 

STL 

240 

19 

65 

MKC 

329 

12 

66 

NEW 

346 

10 

67 

SDF 

MDW 

269 

19 

68 

CAE 

ATL 

194 

18 

69 

BNA 

STL 

254 

17 

70 

MEM 

197 

12 

71 

CLT 

ATL 

227 

14 
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Part  2  -  (Continued) 


Serial 

Number 

Reference 
Community 
Des ignation 

Designation 

Paired  Community 

Distance 

(App.  Statute  Miles) 

Flight 

Exchange/Day 

72 

BHM 

MOB 

1  208 

11 

73 

JAN 

210 

13 

74 

MSP 

OMA 

290 

13 

75 

FAR 

214 

12 

76 

FSD 

198 

12 

77 

LTT 

DAL 

293 

12 

78 

ELP 

TUS 

265 

10 

79 

PHX 

346 

10 

80 

ALB 

BUF 

261 

11 

81 

ABE 

PIT 

239 

11 

82 

FAT 

SFO 

161 

11 

83 

SJC 

138 

10 

84 

JAN 

MEM 

196 

11 

85 

SAV 

ATL 

224 

10 

86 

LAS 

LAX 

234 

93 

87 

ONT 

196 

14 
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Serial 

Number 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 


Part  2  -  (Concluded) 


Reference 

Community 
Des ignation 

Designation 

Paired  Community 

Distance 

(App.  Statute  Miles) 

Flight 

Exchange/Day 

ONT 

SFO 

352 

28 

BO  I 

PDX 

342 

12 

PHF 

RIC 

68 

12 

SYR 

DCA 

300 

14 

SAV 

ATL 

215 

12  , 

SAT 

CVN 

375 

14 

PVD 

DCA 

356 

12 

JAK 

ATL 

285 

30 

ICT 

WDG 

94 

25 

YUM 

FAT 

400 

21 

HFD 

DCA 

304 

10 

Total  exchange  of  flights  per  day 
between  community  pairs  of  Part  2  =  2,477 

Total  exchange  of  flights  per  day 
between  community  pairs  of  Parts  1  and  2  =  5,383 
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Appendix  I 


A  SIMPLE  MODEL  TO  ESTIMATE  THE  FRACTION 
OF  A  ROUTE  FLYABLE  WITH  ROUTE  WIDTH  W 


1 .  Length  of  a  Route  Section  Having  Route  Width  W  As  a  Function  of  the 

Normal  Distance  of  a  Supporting  VORTAC . 

Consider  a  route  AB  and  a  supporting  VORTAC  C  as  shown  in  Figure 
1-1.  Let  the  normal  distance  of  the  VORTAC  from  the  route  by  n  nmi . 

Now  referring  to  Figure  1,  Appendix  D,  page  13  of  the  FAA  Handbook 
7110.18,  ’’Air  Traffic  Control  Service  for  Area  Navigation  Equipped 
Aircraft...",  dated  27  February  1970,  we  have  a  relationship  between 
the  normal  distance  of  a  VORTAC  (referred  to  as  distance  to  point  of 
tangency)  and  the  length  d  of  a  route  section  having  a  width  of  8  nmi 
and  4  nmi.  Some  typical  values  according  to  the  Handbook  are  repeated 
in  Table  1-1  for  convenience.  We  will  use  Table  1-1  to  develop  a  simpli 
fied  expression  giving  the  fraction  of  a  route  flyable  with  route  width 
W  when  the  supporting  VORTAC  grid  is  of  triangular  pattern  with  side 
length.  S . 


A  » 


«  B 


SA- 1096-28 


FIGURE  1-1  A  VORTAC  SUPPORTING  A  ROUTE  SECTION 


Preceding  paSi 
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Table  1-1 


ROUTE  WIDTH  SUMMARY 


Distance  of  Route  Section 
n(nmi) 

Route  Length 
d (nmi ) 

4  nmi 

8  nmi 

0 

50 

100 

25 

44 

96 

50 

16 

88 

75 

— 

70 

100 

— 

20 

2,  Fractions  of  a  Route  Having  Width  W  as  a  Function  of.  the  Side 

Length  S  of  a  Triangular  Grid. 

Consider  Figure  1-2  where  a  triangular  grid  with  side  length  S  is 
shown  along  with  routes  OR,  OQ,  OP,  and  OT  making  angles  of  0,  10,  20, 
and  30°  with  one  of  the  reference  rows  of  the  triangular  grid.  Typically, 
route  sections  having  a  certain  route  width  (e.g.,  4  nmi)  will  consist 
of  small  patches  along  the  route,  as  shown  by  the  shading  in  Figure 
1-2.  The  length  and  distance  between  consecutive  patches  will  be  a 
function  of  the  side  length  S  and  the  angle  between  the  route  and  a 
reference  row  of  the  VORTAC  grid.  The  sum  of  the  lengths  of  patches 
of  a  certain  route  width  will  be  maximum  when  the  route  is  laid  over 
one  of  the  rows  of  the  VORTACs,  e.g.,  along  the  route  OR.  Using  Table 
1-1  as  a  reference,  patches  of  4  nmi  and  8  nmi  were  generated  along 
the  routes  OR,  OQ, . . .  etc.  making  angles  10°,  20°,...  with  S  =  200  and 
S  =  150  nmi  and  assuming  a  route  length  of  1,000  nmi.  The  length  of 
patches  of  4  nmi  and  8  nmi  route  widths  were  summed  along  each  route 
and  divided  by  the  total  route  length,  i.e.,  1,000  nmi,  so  that  along 
each  route  the  totals  of  the  sections  having  a  width  of  4  nmi  and  8  nmi 


FIGURE  1-2  ROUTE  SECTIONS  WITH  CERTAIN  ROUTE  WIDTHS 


can  be  expressed  as  fractions  of  the  route  length.  Some  typical  values 
are  shown  in  Table  1-2. 

Table  1-2 

FRACTIONS  OF  A  ROUTE  WITH  WIDTHS  OF  4  NMI  AND  8  NMI 


Route  Angle  with 
Reference  Row 

S  =  200  nmi 

S  -  150  nmi 

4  nmi 

8  nmi 

4  nmi 

8  nmi 

0 

.25 

.50 

.35 

.65 

±10 

.  20 

.44 

.30 

.62 

±20 

.18 

.42 

.28 

.58 

±30 

.115 

.36 

.20 

.52 

±40 

.18 

.42 

.28 

.58 

±50 

to 

o 

.44 

.30 

.62 

±60 

.25 

.50 

.35 

.65 

139 


The  relationships  in  Table  1-2  are  depicted  graphically  in  Figure 
1-3.  Note  that  the  fractions  for  routes  making  10° and  50°  angle  with 
a  reference  row  are  identical,  since  a  route  that  is  inclined  at  50° 
to  one  row  is  actually  inclined  at  10°  to  another  row  of  the  grid. 
Similarly,  20°  inclination  is  equivalent  to  a  40°  inclination  and  so 
on.  This  is  because  of  the  repetitive  nature  of  a  triangular  VORTAC 
grid  (every  60°). 


z 

o 


FIGURE  1-3  ROUTE  WIDTH  SUMMARY  FOR  ROUTE  SECTIONS  (S  =  150  AND  S  =  200  nmi) 


Without  any  significant  error,  a  linear  fit  can  be  made  for  the  data 
points  as  shown  in  Figure  1-3.  If  we  now  assume  that  for  other  values 
of  S  and  W,  similar  fits  can  be  established  using  linear  relationships, 
i.e,,  if  side  length  S  is  increased  by  a  factor  k^,  the  fractions 
for  a  certain  route  width  will  decrease  by  a  factor  k^--note  that 
larger  values  of  S  mean  a  smaller  number  of  patches  of  a  certain  route 
width  for  a  fixed  route  length — and  if  W  is  chosen  higher  or  lower  a 
corresponding  increase  or  reduction  in  the  fractions  would  occur,  then 
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using  the  S=200,  W  =  4  nmi  set  of  data,  a  relationship  giving  the  fraction 
of  a  route  having  a  width  W,  side  length  S,  and  route  angle  0  with  one 
of  the  reference  rows  of  the  grid  can  be  expressed  analytically  as 
follows : 


F 

V 


200 

s 


0.25 


.  25 (W  -  4) 
4 


125  e 
30 


for  0°  £  l9l  ^30° 
and  60°  ^  1 6 1  ^  90° 


200 

S 


0.25 


,  25( W  -  4) 
4 


.125(60  -  9) 
30 


for  30°  ^  lO!  ^  60° 


where  F  =  fraction  of  the  route  having  route  width  W. 

V 

The  above  relationship  is  of  course  not  exact,  but  seems  to  be 
adequately  correct  for  a  gross  evaluation  of  matching  between  a  route 
network  and  a  given  triangular  V0RTAC  grid. 
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